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ABSTRACT 


Stress  fields  in  the  components  of  a  multifiber,  unidirectional  com¬ 
posite  that  are  caused  by  shear  loading  in  axial  directions  are  studied  and 
solutions  are  presented.  Computations  for  variou?  volumetric  contents  and 
material  constants  are  described,  as  well  as  experiments  which  showed  good 
correlation  with  the  theoretical  results  from  this  work.  Further,  it  is 
shown  that  certain  mathematical  solutions  normally  used  for  such  problems 
are  not  compatible  with  present  boundary  conditions,  and  that  the  problem 
becomes  two-dimensional. 
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INTRODUCTION 


The  problem  of  shear  loading,  together  with  axial  and  transverse  loading, 
constitutes  the  general  case  of  oblique  loading  of  a  composite,  A  previous 
contract,  DA  44 -177-AMC -44 1 (T) resulted  in  the  solution  of  the  axial  load¬ 
ing  case  of  a  finite-length,  multifiber  composite.  The  model  selected  was 
a  nuitifiber  composite  of  hexagonal  arrangement,  loaded  in  shear  in  such  a 
manner  as  to  represent  an  external  attack  on  surfaces  parallel  to  the  fibers. 

The  computer  program  obtained  takes  into  account  the  fiber  and  matrix 
modulus,  the  fiber  radius,  and  the  fiber  volumetric  content.  To  obtain  the 
stresses,  analytical  solutions  have  been  used  and  adjusted  to  the  boundary 
conditions  at  the  interfaces  between  fibers  and  matrix,  to  certain  repeating 
symmetry  conditions,  and  to  otner  geometrical  conditions  characterizing  the 
problem.  The  boundary  conditions  were  matched  in  368  points  for  each  fiber 
and  the  hexagonally  surrounding  resin. 

The  composite  longitudinal  shear  modulus  obtained  for  a  fibrous  com¬ 
posite  was  compared  with  the  results  from  previous  work  done  in  this 
field.  A  design  formula  has  been  derived. 
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TECHNICAL  DISCUSSION 


The  geometrical  arrangement  of  the  reinforcements  in  the  matrix  is 
assumed  to  be  such  that  a  central  fiber  is  surrounded  by  six  symmetrically 
spread  fibers  which,  in  turn,  are  surrounded  by  six  others.  Figure  1  shows 
the  plane  geometry  of  the  composite. 


r 


Figure  1.  Plane  Geometry  of  the  Composite  With 
the  Characteristic  Element  Shaded. 


Figure  2  shows  a  composite  before  and  after  shear  deformation. 

It  is  assumed  that  the  material  constants,  such  as  modulus  of  elasticity 
and  Poisson's  ratio  of  both  the  reinforcement  and  the  matrix,  shall  differ 
from  one  another.  The  transmission  of  loads'  from  one  material  into  the 
other  is  expressed  by  the  boundary  condition  equations  at  the  interfaces. 

The  kinematic  boundary  conditions  will  be  assumed  as  complete  contact  be¬ 
tween  the  two  materials  so  that  no  slippage  will  be  possible. 
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Figure  2.  Composite  Before  and  After  Deforming. 

As  in  the  axial  loading  case,*  an  attempt  was  made  to  express  the  solu¬ 
tions  of  the  Navier  equations  in  terms  of  Neuber -Papkovitch  potentials 

.  These  potentials  are  defined  as  being  solutions  of  the  Laplace  equa¬ 
tion 

^P.  =  0  (1) 

where  i  =  0 , 1 , 2 ,3  . 

In  equation  (1),  P0  is  scalar  while  the  other  terms  are  vector  com¬ 
ponents.  In  rectangular  coordinates,  as  shown  in  Figure  3,  the  displace¬ 
ment  solutions  of  the  Navier  equation  are 


nrw  sr  (xJpj +  p°> 


where  i 

j 


x,y,z  and 
x,y,z  . 


*  Contract  DA  44 -l 77-AMC-44 1 (T) . 
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Figure  3.  Constant  From  the  Hexagonal  System  for  Definition 
of  the  Coordinates  Used  in  the  Analysis  and  in 
the  Computation. 


Applying  a  shear  X  on  the  surface 
following  symmetry  conditions  would 

x  =  ±A  ,  one  would  assume  that 
be  adequate: 

the 

ux  (x,y,z) 

*  "ux  <-x>-y>z) 

(3) 

ux  (x,y,z) 

=  ux  (x,y,-z) 

(4) 

'4 


(x.y.z) 

=  -uy  (-x , -y ,z) 

(5) 

“y 

(x,y,z) 

=  (x,y,-z) 

(6) 

uz 

(y  ,y  »z) 

=  ~u2  (-x,-y,z) 

(7) 

(x ,y,z) 

=  uz  (x,y, -z) 

(8) 

But  because  the  solutions  of  equation  (1)  are  of  the  form 


sin  ,  .  sin  .  sinh  ,  N 

(ax)  *  (By)  *  ,  (yz) 

cos  cos  cosh 


it  is  impossible  to  find  in  equation  (9)  a  combination  of  solutions  which 
satisfies  the  symmetry  conditions  represented  by  equations  (3)  through  (8). 
On  the  basis  of  physical  considerations  alone,  there  can  be  no  objection 
to  these  symmetry  conditions.  The  problem,  however,  stems  from  the  fact 
that,  for  large  displacements,  the  surfaces  x  =  ±A  would  try  to  come 
closer  to  each  other  during  the  application  of  the  shear  T  .  Since,  in 
the  theory  of  elasticity,  only  infinitesimally  small  displacements  are 
taken  into  account,  the  surfaces  x  =  ±A  do  not  come  closer  while  shear 
is  applied.  The  interpretation  of  this  apparent  contradiction  is  that  it 
is  impossible  to  find  an  equilibrated  and  compatible  solution  which 
satisfies  the  symmetry  conditions.  In  reference  to  the  present  problem, 
it  can  be  stated  that  the  solutions  of  equation  (9)  will  never  satisfy  the 
symmetry  conditions  in  equations  (3)  through  (8). 

In  the  following  discussion,  symmetry  conditions  will  be  presented 
which  comply  with  the  solutions  of  equation  (9).  Nevertheless,  a  pre¬ 
requisite  to  further  use  of  these  symmetry  conditions  (which  will  now  be 
compatible  with  the  solutions  of  elasticity)  is  the  determination  of 
whether  or  not  the  proposed  symmetry  conditions  are  adequate  for  the 
problem.  A  set  of  symmetry  conditions  compatible  with  the  solutions  is 
given  below. 


[ 

ux  (*,y,z) 

=  ux  (-x,-y,z) 

(10) 

ux  ( x,y,z ) 

=  -ux  (x , y , -z) 

(ID 

u  (x,y,z) 

=  u  (-x,-y,z) 

(12) 
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uy  (x,y,z)  =  -Uy  (x,y, 

-z) 

(13) 

uz  (x,y,z)  =  -uz  (-x,- 

y.z) 

(14) 

u  (x,y,z)  =  u  (x , y , - 

z  z 

z) 

(15) 

Proceeding  then  with  conventional  procedures,  we  must  select  such  combina¬ 
tions  of  equation  (9)  which,  when  introduced  into  equation  (2), satisfy  the 
conditions  in  equations  (10)  through  (15).  These  combinations  are  presented 
below. 

=  (Aj  sin  ax  sin  By  +  Ag  cos  ax 

cos  By)  sinh  \z 

(16) 

P  =  (Bi  sin  ax  sin  By  +  Bg  cos  ax 

y 

cos  By)  sinh  yz 

(17) 

P  =  (Cx  sin  ax  cos  By  +  Ca  cos  ax 

z 

sin  By)  cosh  yz 

(18) 

P0  =  constant 

In  equations  (16)  through  (18),  A  ,  B  ,  C  ,  a  ,  and  6  are  independent  con¬ 
stants,  while  y  =  aS  +  @2  .  Additionally,  from  the  geometry  shown  in 
Figures  1  and  3,  it  could  be  assumed  that,  at  the  surface  x  =  +A  ,  the 
displacements  uz  in  axial  direction  z  are  constant. 

Also  constant  at  this  surface  is  the  displacement  in  the  x  direction. 
The  shear  in  the  y  direction  in  the  plane  x  =  ±A  is  zero.  Consequently, 
we  have  the  following  boundary  conditions  at  the  planes  x  =  ±A  : 

uz  (A,y,z)  = 

Ki 

(19) 

uz  (-A,y,z)  = 

-Ki 

(20) 

c 

X 

> 

N 

tw/ 

II 

Ka 

(21) 

ux  (-A’y’z)  = 

-Kb 

(22) 

CT*y  (A)  =  axy  (‘A)  = 

0 

(23) 
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Using  Figure  3,  we  can  visualize  that,  at  the  surface  y  ■  ±B  , 


a  (x,B,z)  ■  0 
y* 

(24) 

(x ,  - B ,  z )  =  o 

(25) 

a  (x,B,z)  =  0 

y4 

(26) 

CT  (x,-B,z)  =  0 

yz 

(27) 

uy  (x,B,z)  =  K3 

(28) 

uy  (x  ,  -  B,  z)  =  -Ka 

(29) 

antisymmetric  functions  in  z 

(equations  11  and  13) 

Since  u  and  u  a 
x  y 

then  Ka=0  ,  on  the  basis  of  equations  (21.)  and  (22),  and  ^*=0  ,  from 

equations  (28)  and  (29).  Consequently,  the  only  non homogeneous  boundary 

condition  is  u  =Kj  .  If  only  integral  boundary  conditions  are  imposed  at 

both  free  ends,  then  these  bourJaries  do  not  contribute  to  the  solution  of 

the  problem;  they  are  automatically  satisfied  Because  of  the  equilibrium 

conditions  and  because  there  are  no  applied  fu'-ces  of  known  distribution  at 

that  surface.  Since  the  fiber's  length  is  large  compared  with  its  diameter, 

it  is  permissible  to  assume  that  the  displacements  u  ,  u  ,  and  u  do 

x  y  z 

not  depend  on  z  .  This  assumption  is  valid  for  the  entire  fiber  and 
becomes  invalid  only  on  one  or  two  fiber  diameters  from  the  fiber  ends. 

(For  the  ends,  a  perturbation  analysis  will  be  performed  later.)  With  this 
assumption,  the  displacements  in  equation  (2)  will  be  reduced  to  the 
following  for-": 


1  * 

UX  ■  PX  •  40^0  *7  <xF*  +  yPy) 

uy  ■  F»  -  ^  (xP» +  yV 

u  =  P 


(30) 

(31) 

(32) 
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From  the  displacements  in  equation  (2)  or,  more  specifically,  in  equations 
(30)  through  (32),  the  strains  are  obtained  by  the  relation 


'ij 


rUl ,  3u, 

— L  + 

*x.  dx 
3-  j 


(33) 


Hook's  law1  furnishes  the  relation  between  stress  and  strain: 


'ij  1+v  ^£ij  +  6ij  l-2v  SkkJ 


(34) 


From  equations  (33)  and  (34),  we  obtain 


ij  2  (1+V) 


2v 

dx.  cix,  ij  l-2v  dx 

l  j  J  k 


(35) 


In  order  to  calculate  the  stresses  in  terms  of  Neuber-Papkovitch  potentials, 
the  expressions  for  the  displacements  given  in  equation  (2)  will  be  intro¬ 
duced  into  equation  (35),  as  follows: 


ij  2  (1+V) 


4 (1  -v)  (xkPk  +  Po>}  + 


_a 

dx 


[Sxi  b 

;{Pi  -  4(1^)  ^  <*<  Pt  +  P=  j 


6lj  1^7  ^  {P-  -77T^7^-(-.p.+  ro) 

J  m 


m  4(l-v)  ^x  n  n 
m  *  m 


1] 


(36) 


or 


2  (1+v) 


*P.  ?P, 


+  6 


2v  m 


3x.  5x ,  ij  l-2v  ^x  2(l-v)  Sx.ox.  k  k 

i  J  m  i  J 


(*,A  +  P0  )  - 


2v 


ii  4(l-2v)(l-v)  *x  n  n 

m  m 


(x  P  +  P0  ) 


(37) 
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By  using  engineering  notations  a  ,  cr  ,  o  >  and  so  forth,  we  can  write 
equation  (37)  ii.  another  form:  y  y 


a  =  ct 

XX  X 


_  I  U  i  .  I  * 

E  _ x  1  _ x  _ 

2 ( 1+v)  l  Ax  ‘  2C1-V)  X  ^  y  ^ 


i )  +  _y_  !Li] 
;  I  1-v  ay  J 


E 

r  ap  i 

y  i 

A2  P  a2  P 

x  +  y  ^ 

AP  1 

+  _ 21 

°yy  "  °y  2(1+V) 

Ay  ‘  2 ( 1  -V ) 

■  Ay2  Ay5 

+  1-V  Ax  J 

2  (1+v)  (1 


v _ 

-v)  . 


*PX 

*x  +  *  \ 


l-2v 

/AP 

AP  \ 

—  + 

__Z 

2(1 -V) 

» 

*y 

Ax  J 

2(1 -V) 


[X  AxAy  +  ^  AxAy 


E  z 

2  ( L+v)  Ax 


E  z 
2  (1+v)  Ay 


Some  important  conclusions  characterizing  tno  shear  problem  can  be  made 
from  the  stress  equations,  (38)  through  (43)  In  these  equations,  (38) 
through  (41)  depend  only  on  P  and  P  ,  and  (42)  and  (43)  depend  only 

on  P According  to  the  theorem  that  a  harmonic  function  is  zero  if 
the  boundary  conditions  are  completely  homogeneous,  the  potentials  P 
and  P  introduced  into  the  homogeneous  boundary  conditions  expressed 

in  equations  (21)  through  (27),  keeping  in  mind  that  K2  is  also  zero, 
must  be  zero.  Therefore, 

P  (x ,  y ,  z)  =  P  (x,y,z)  =  0  (‘ 

x  y 
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Consequently,  during  shear  application  the  following  displacements  and 
stresses  are  zero,  not  only  at  the  boundaries  but  also  in  the  whole  com¬ 
posite  : 


(45) 


In  conclusion,  this  problem  must  be  recognized  as  two-dimensional ,  having 
only  the  following  displacement  and  stresses: 


u 

z 


a 


xz 


and  a 

yz 


Remembering  the  definition  of  the  potentials  P  (equation  1)  and  their 
relationship  with  the  displacements  (equation  2)  and  the  stresses  (equa¬ 
tions  42  and  43),  the  longitudinal  shear  problem  can  then  be  reduced  to 
the  solutions  of  the  Laplace  equation  and  their  adjustment  to  the  boundary 
conditions . 


=  0 
z 


Furthermore,  since  uz  =  pz  » 


V2  u  =  0  (46) 

z 

The  expressions  for  the  shear  stresses  are  then  expressed  by  the  well-known 
relation 


a 


xz 


E 

2 (1+v) 


3u 

_ z 

3x 


(47) 


o 

yz 


E 

2 (1+v) 


1^8) 


The  displacements  and  stresses  must  satisfy  certain  symmetry  and  boundary 
conditions,  which  are  restated  in  the  following  paragraphs. 
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SYMMETRY  CONDITION 


BOUNDARY  CONDITIONS 


u  (x,y)  -  -u  (-x,-y) 

z  * 


GyZ  (X»B>Z)  "  0 


CTyz  (x>"B>z)  =  0 


uz  (A, y , z)  =  Ki 


uz  (-A,y,z)  =  -Kx 


Along  the  interfaces,  the  contact  shear  stresses  must  be  continuous  in 
both  materials  I  and  II  (fiber  and  resin). 


cr^  cos  cp  +  sin  cp 
xz  T  yz 


II  II 

a  cos  cp  +  ct  sin  cp 
xz  T  yz 


Expressed  by  equations  (47)  and  (48),  the  boundary  condition  in  equation 
(54)  is 


du  du  TT  9u  ,,  du 

G1  ~  cos  cp  +  G  -jy  sin  cp  =  G11  —  cos  9  +  G  sin  cp 


where 


2  (1+V  1,11  ) 


and  where  the  superscripts  I  and  II  designate  reinforcement  and  matrix 
respectively. 
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The  first  approach  tried  was  the  finite  difference  method,  using  a 
square  net  of  112  nodes,  a  Laplace  operator  of  second-order  approximation 
[0  (h  )]  ,  and  a  linear  interpolation  along  the  interface  curved  boundary. 

The  following  problem  appeared  in  the  solution  obtained  for  a  test  case 
with  65%  reinforcement. 

In  the  fiber,  near  the  interface,  the  displacement  function  in  the 
z  direction  u  had  a  convex  curvature  with  respect  to  the  displacement 
in  t;ie  resin.  This  is  illogical,  for  the  following  reason;  Close  to  the 
surface,  where  the  transition  takes  place  from  one  material  into  the  other, 
the  stress  must  increase  rather  than  decrease,  as  a  convex  curvature  would 
indicate.  This  is  shown,  in  exaggerated  form,  in  Figure  4. 


Figure  4.  Shear  Displacement  in  the  Axial  Direction, 
Obtained  by  the  Finite  Difference  Method 
and  Compared  With  the  Point -Matching  Method. 
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vff'  ■ 


Since  the  modulus  of  elasticity  is  calculated  from  the  relaLion 


a 

xz 


(56) 


that  will  say  from  the  slope  of  the  displacement  uz  as  shown  later  in 
equation  61,  the  two  methods  will  give  two  different  moduli  of  elasticity, 
corresponding  to  the  slopes  of  the  two  curves  in  Figure  4 ,  in  x  =  A.  The 
relation  between  the  slope  of  the  displacement  in  the  reinforcement  and 
that  in  the  resin  is  obtained  from  the  fact  that,  at  the  interface,  the 
shear  stresses  are  equal  in  both  materials: 


G 


I 


or 


(57) 


(58) 


This  relation  can  be  very  large  in  number.  In  composite  materials,  it 
ranges  from  10  to  180;  it  is  therefore  difficult  to  represent,  let  alone 
satisfy,  the  boundary  conditions  at  the  interface  with  the  finite  differ¬ 
ence  method.  In  the  point-matching  method,  different  functions  for  both 
materials  will  be  used. 

METHOD  OF  SOLUTION 


A  boundary  collocation  (point  matching)  method  is  used.  This  method 
essentially  consists  of  choosing  a  set  of  linearly  independent  functions 
which  satisfy  the  differential  equation  exactly.  The  total  solution  com¬ 
prises  the  solution  of  individual  functions,  each  one  multiplied  by  a 
free  constant.  These  constants  are  determined  by  postulating  that  the 
boundary  conditions  at  a  certain  number  of  points  must  be  satisfied.  If 
the  number  of  boundary  conditions  is  equal  to  the  number  of  undetermined 
coefficients  in  the  functions,  a  system  of  linear  equations  with  exactly 
as  many  equations  as  unknowns  will  result.  If,  however,  more  boundary 
points  are  used  in  the  equations  than  there  are  free  constants,  an  "over¬ 
determined"  system  of  linear  equations  will  arise,  which  can  be  solved  by 
the  classical  least-squares  method. 
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With  boundary  collocation  methods,  it  is  usually  safer  to  specify  more 
boundary  conditions  (about  twice  as  many  seems  to  be  a  good  rule  of  thumb). 
This  results  in  a  smoother  overall  solution. 

The  success  of  any  collocation  method  depends,  to  a  large  extent,  on 
a  suitable  choice  of  the  set  of  functions  by  which  the  solution  is  to  be 
approximated.  The  following  functions  would  satisfy  both  the  Laplace 
equation  and  the  symmetry  condition: 


and  a  similar  function  for  w  (all  superscripts  in  the  function  are  re¬ 
placed  by  II). 


These  functions  were  programmed  first,  but  their  oscillatory  nature  did 
not  allow  a  reasonable  fit  of  the  boundary  conditions  along  the  line  x  =  A  . 


A  second  set  of  functions,  which  also  satisfies  the  differential  equa 
tions,  the  symmetry  conditions,  and,  in  addition,  the  boundary  conditions 
along  y  ■  ±B  ,  is 


u2(x,y) 


e 


m*l,2,3L 


.  I  .  ,  m  tt 
A  sinh  —r—  x 
m  B 


Im  tt  ,  „I  ,  / (2m-l )  tt 

cos  -  y  +  B  cosh  P — - — i — 

B  3  m  \  2B 


x  X 


I  (2m-l )  tt 

sin  2B  y 


•+  E^x 


and  a  similar  function  for 


II 

w 


These  functions  have  the  same  basic  oscillatory  character  and  are 
therefore  not  suitable  as  a  solution  to  the  shear  problem. 

The  functions  which  were  finally  chosen  are  polynomials,  and  u^(x,y) 

II  2 

and  u^  (s,y)  were  assumed  to  be  of  the  following  form.  These  functions 

are  mentioned  because  they  seem  to  be,  but  are  not,  usable  for  the  collo¬ 
cation  method.  The  present  case  is  one  of  the  examples  which  has  shown 
chat  not  all  the  functions  that  satisfy  the  Laplace  equation  and  the  symmetry 
conditions  can  be  used  successfully  in  a  collocation  method. 
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u‘(x,y) 


max 

E 

j-i 


aj  Re  (x  +  iy)(2j_1^  +  b*  Im(x+ty) (2 j _1 ^  + 


max 

E 

k  =  1 


c}  Re  (x  +  iy)2^  +  Im(x  +  iy^2‘C^  +  eJ 
K  K 


(59) 


II,  ^ 
uz  (x,y)  = 


max 


Y*  a11  Re  (x  +  iy)(2m_1)  +  bH  I  (x  +  iy)2m_1 

4"  "  "  ‘ 

m  =  1 


m  m 


(60)* 


The  above  expression  for  u  does  not  satisfy  the  symmetry  condition. 

Since,  in  the  basic  element,  the  two  areas  representing  fibers  are  not 

connected,  it  is  not  necessary  that  the  same  function  be  used  in  the  two 

fiber  areas.  The  function  u^*  ,  which  extends  over  a  connected  area  both 

2 

to  the  right  and  to  the  left  of  the  y-axis  and  above  and  below  the  x-axis, 
must  satisfy  the  symmetry  condition  ull(x,y)  =  -u^  (-x , -y)  .  The  above 

function  has  this  property. 

The  appendix  gives  details  on  the  boundary  points  where  the  boundary 
conditions  were  established. 


*  This  function  was  proposed  by  Dr.  George  Burgin  of  Decision  Sciences. 
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f  Displacements  When 
s  Applied  at  a  Sur- 
to  the  Surface  +A. 


aceraents  in  a  composite  d 
+1  displacement  at  the 


1 ac  ement 
Boundary 


STRESS  DISTRIBUTIONS 


Figures  9  through  12  show  the  stress  distribution  in  the  composite  com¬ 
ponent  that  is  caused  by  longitudinal  shear  at  the  surface  x  =  ±A  .  In 
these  representations,  the  stresses  generated  are  divided  by  the  applied 
stress,  so  that  the  numbers  shown  are  the  shear  stress  concentrations.  The 
applied  stress  is  expressed  by  the  average  applied  load  at  the  boundary 
surface 


cp,  deg 


Figure  9.  Shear  Stress  at  the  Interface  due  to 

Longitudinal  Applied  Shear  Load  axz(A)* 
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60%  Boron 


0.532  0.799  1.065  1.352  1.599 


Shear  Stress  aXz 


as  Function 


y 


23 


■s**?'-"***- 


LONGITUDINAL  SHEAR  MODULUS  OF  A  COMPOSITE 

The  total  shear  modulus  of  the  composite  can  be  calculated  from  the 
applied  shear  Load  necessary  to  produce  the  shearing  angle  y  .  Expressing 
the  total  she  ir  load  T  at  the  surfaces  x  =  ±A  by  the  shear  stress 
produced  at  the  surfaces  in  both  materials,  we  obtain 


■r 


(a  )  idy 
xz 

x=A 


J  (Oxz)n  4<ly  +  J  (^xz)1  idy 
4*  x=A  I da  X"A 


B-a  /A u  \  •Li  -B  ib u  V1 

/  kL tdy  +  G  /  ''dy 


In  addition  to  the  expression  in  equation  (62),  the  total  shearing  force 
of  the  composite  can  be  expressed  by  the  composite  shear  modulus  G: 


T  6 

2tB  '  ‘  G  I 


In  a  hexagonal  arrangement, 


A  =  B 


■F~ 


,  i  /7  i- 

2  V  tl 


Consequently,  the  composite  shear  modulus  is 


■  *4-  r  « *  "■  /  et.  •> 

^B  B-a 
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From  the  computer  program,  as  described  in  the  Appendix,  the  composite 
shear  modulus  G  has  been  presented  as  the  function  of  the  fiber  modulus 
Gf,  the  matrix  modulus  Gm,  and  the  volume  percentage  of  the  fiber  Vf 
(see  Figure  13). 


0.72 

0.70 

0.65 

0.60 

0.55 


Figure  13.  The  Longitudinal  Shear  Modulus  of  ;i 
Unidirectional  Composite  as  Computed 
From  the  Computer  Program. 


In  the  following  test,  the  computer  results  are  compared  with  the  results 
of  the  approximation  formula  by  Rosen,  Foye ,  Ekvall  and  Berg, 2-5  the  com¬ 
puter  results  of  Tsai,1'  and  the  test  results  of  Adams  and  Doner7.  The 
results  calculated  from  Rosen's  formula*" 


G_ 

Gm 


jr-  (l  +  vf )  +  (i  -  vf ) 

m _ 

Gf 

~  (i  -  vf )  +  (l  +  vf ) 

m 
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a8ree  <luite  "e11  °this°oneCis  SmS  compared  with  our  oom- 

;it «“  £— ua 


_G_ 

Gm 


are  287.-347.  lower  than  our  computer  results 


O  Rosen's  Formula 
—  Computer  Results 


Gf/Gm 


Figure 


14.  Computer  Results 


Compared  With  Formula  by  Rosen 


The  results  from  the  equation  given  by  Ekval 


G/G 


Figure  16  is  a  comparison  of  our  computer  results  with  the  results 
of  Foye's  formula.3 


8.0 


7.0 


Computer  Results 
Formula  by  Foye 


Figure  16.  Computer  Results  Compared  With  Formula  by  Foye. 
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Test  results  obtained  by  Adams  and  Doner  are  slightly  higner  but  very 
close  to  the  computer  results  obtained  under  this  program  (Figure  17). 

An  approximate  equation  deveijped  from  the  present  computer  program  is 


~  =  (11.78  V  2  -  13  V  +  3.78)  tnl^r  6.66  V  -  1.27  (66) 

Cm  f  f  \Gj  f 


Figure  18  shows  equation  (66)  in  "omparison  with  the  computer  results 
obtained  in  this  work. 


EXPERIMENTS 


The  experimental  portion  of  this  contract  was  conducted  primarily  to 
provide  clarification  of  the  stress  concentrations  at  the  'id  of  the  fibers. 

For  this  purpose,  two  test  specimens  such  as  those  shown  in  Figure  19 
were  fabricated.  The  resin  used  was  epoxy  type.  Glass  fibers  measuring 
0.28  inch  in  diameter  were  embedded,  in  a  hexagonal  array,  into  the  resin. 
Using  an  epoxy  resin  as  the  adhesive,  the  specimen  was  attached  to  a  steel 
test  frame,  as  illustrated  in  Figure  20.  With  this  frame,  shear  loads 
could  be  applied  to  the  specimen.  The  specimen  was  analyzed  by  photo¬ 
elastic  methods,  first  without  load  and  then  with  incrementally  increasing 
loads.  Figures  21  and  22  show  the  specimen  illuminated  by  polarized  light, 
in  the  unloaded  state  and  loaded  with  800  pounds  of  pressure,  respectively. 
It  is  possible  to  see  only  small  stress  concentrations  at  the  end  of  the 
fibers.  The  depth  of  this  change  in  the  stress  is,  more  or  less,  two 
fiber  diameters  along  the  axis. 

In  conclusion,  the  stress  analysis  with  its  assumptions  made  for  long 
fibers,  as  developed  during  this  contract,  is  justified.  The  results  are 
exact  along  the  fiber,  with  the  exception  of  a  very  small  part  close  to 
the  end.  In  practical  use,  the  length  of  this  region  (no  more  than  four 
fiber  diameters)  is  about  1/2000  for  the  shortest  fiber.  Therefore,  the 
influence  of  stress  concentrations  at  the  fiber  end  is  negligible  for  the 
computation  of  composite  shear  modulus. 
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Section 
through 
the  sped 


H 


Test  Picture  Frame  for  the 
Shear  Te st  . 
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CONCLUSIONS  AND  RECOMMENDATIONS 


It  appears,  from  this  analysis  of  a  composite  material  subjected  to 
longitudinal  shear  loading,  that  micromechanical  stresses  occur  which  are 
larger  than  the  stresses  applied  to  the  composite  as  a  whole.  This  report 
illustrates,  for  example,  that  the  stress  Cnz  along  the  interface  is  357. 
larger  than  the  applied  average  stress  axz(A)  for  an  angle  $  =  28° 
(reference  Figure  9).  In  other  words,  a  stress  concentration  exists  which 
is  induced  by  the  inclusion  of  fibers  in  the  matrix.  This  form  of  stress 
concentration  can  be  known  only  by  using  the  micromechanical  analysis 
based  on  the  theory  of  elasticity,  as  used  in  this  report. 

For  the  most  part,  the  shear  modulus  of  elasticity  obtained  by  this 
micromechanical  analysis  correlated  well  with  test  results.  It  must  be 
noted,  however,  that  obtaining  the  shear  modulus  by  experimental  methods 
(torsioned  tubes,  plates,  etc.)  is  difficult,  from  the  standpoint  of 
accuracy.  This  was  the  cause  of  the  variation  in  the  results  reported. 

The  formula  for  shear  modulus  which  emerged  from  this  work  yields  a 
high  degree  of  accuracy  for  most  composites  used  in  structures.  This 
formula  can  be  used  easily  by  the  designer. 

Once  the  internal  behavior  of  composite  materials  under  load  is 
understood,  then  it  will  be  possible  to  attain  material  optimization  by 
judiciously  selecting  the  most  adequate  combination  of  geometry  and  mech¬ 
anical  properties  of  the  fibers  and  the  matrix.  This  contract  effort  was 
directed  to  this  end.  Work  is  being  continued  on  the  micromechanics  pro¬ 
gram  in  the  area  of  other  forms  of  loading,  such  as  transverse  loading. 
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APPENDIX 


NUMERICAL  CALCULATIONS  AND  COMPUTER  PROGRAM 


The  program  presently  uses  the  following  upper  summation  limits: 


JMAX  =  KMAX  =  MMAX  =  10 


This  gives  a  total  number  of  61  unknown  coefficients.  Ninety-one  condi¬ 
tions  at  the  boundary  points  are  specified,  so  that  a  linear  system  of 
91  equations  with  61  unknowns  must  be  solved. 

The  above  expression  for  u^  equation  (59)  does  not  satisfy  the  sym¬ 
metry  condition.  Since  in  the  basic  element  the  two  areas  representing 
fibers  are  not  connected,  it  is  not  necessary  that  the  same  function  be 
used  in  the  two  fiber  areas.  The  function  w*I  ,  which  extends  over  a 
connected  area  both  to  the  right  and  to  the  left  of  the  y-axis  and  above 
and  below  the  x-axis,  has  to  satisfy  the  symmetry  condition  wll(x,y)  = 
-wII(-x ,-y)  .  The  function  given  in  equation  (60)  has  this  property. 

Figures  23  and  24  show  the  points  u;  '  r  the  \ oundary  collocation 
for  the  two  different  cases  a<A  and  a>A  .  Tables  I  and  It  show  the 
91  equations  generated  for  both  cases. 

DESCRIPTION  OF  THE  PROGRAM 


The  program  is  written  in  FORTRAN  63  (similar  to  FORTAN  IV)  for  a  Con¬ 
trol  Data  Computer  Model  3600.  It  consists  of  approximately  1200  FORTRAN 
statements.  The  following  general  flow  chart  gives  an  idea  about  the 
options  available  in  the  program. 

THE  SUBROUTINES  AND  THEIR  PURPOSES 

ALBERTO  This  is  the  main  program.  It  reads  and  prints  the 
input  data.  It  calls  the  necessary  subroutines  to 
calculate  the  equations  for  the  coefficients,  to 
solve  these  equations,  to  calculate  displacements 
and  stresses;  and  to  calculate  the  ratio  of  the  two 
shear  moduli 

Slatr  ix/^Total 

GETMTRX  This  routine  generates  equations  at  91  boundary 
points.  Uses  subroutine  power. 
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Interface 


Figure  23.  Arrangement  of  Boundary  Points  for  a  <  A 


TABLE  I  EQUATIONS  GENERATED  IF  p<a 


Point 


Cond i t ion 

No .  o  £ 
Poin  ts 

No .  of 
Equation 

I 

11 

z 

II 

=  u 

z 

1 

] 

I 

T 

xz 

II 

=  T 

xz 

1 

2 

_I 

yz 

=  0 

1 

3 

II 

T 

yz 

=  0 

1 

4 

T 

U  * 

Z 

II 

-  u 

z 

14 

5-18 

T 

T 

XZ 

COScp  + 

^  yz 

sincp  = 

14 

19-32 

II 

T 

XZ 

.  II 

cosm  +  t 

yz 

sincp 

I 

u 

z 

=  k 

1 

33 

II 

u 

z 

=  k 

1 

34 

I 

T 

yz 

II 

yz 

1 

35 

yz 

=  0 

9 

36-44 

I 

'Y 

yz 

=  0 

1 

45 

I 

u 

z 

=  k 

1 

46 

TABLE  I  Continued 

Point 

Condition 

No .  of 

Points 

No.  of 
Equat i on 

f 

u1-  k 
z 

9 

47-55 

g 

11 

u  =  k 
z 

4 

56-5^ 

h 

11  =  1C 

1 

60 

h 

T11  =  0 
yz 

1 

61 

i 

11  n 

V"  0 

12 

62-/3 

j 

Tn=  0 
yz 

3 

74-76 

f 

QuI/9y  =  0 

2 

9 

77-85 

c 

auI/5y  =  0 
z 

1 

86 

c 

du^/dy  =  0 

1 

87 

g 

du^/dy  =  0 

4 

88-91 

42 


TABLE  II.  EQUATIONS  GENERATED  IF  p>a 


Po  in  t 

Condition 

No .  of 

Points 

No .  of 
Equat ion 

<J 

T1  =  0 
yz 

12 

36-47 

e 

T1  -  0 
yz 

1 

48 

e 

u1-  k 
z 

1 

49 

f 

c 

II 

9 

50-58 

g 

11  L. 

u  =  k 

z 

4 

59-62 

h 

11  L. 

t  =  k 

z 

1 

63 

h 

11  n 
yz 

1 

64 

i 

n  _ 

T  =0 

yz 

12 

65-76 

NOTES:  1,  Equations  1  through  35  are  the 
same  as  for  the  case  a  <A  . 


2.  There  are  no  points  j 


3.  Equations  77  through  91  are  the 
same  as  for  the  case  a  <  A. 


POWER 


AM  IN  US 


PWRF 


NORMEQ 


LISTARAY 


riMATINV 


This  subroutine  computes  both  the  real  and  the 
imaginary  parts  of  the'  function 

(x  +  iy)n  ,  where  n  S20 

The  routine  is  first  entered  through  the  entry 
point,  initial,  in  order  to  calculate  binomial 
coefficients.  The  entry  points  POWER,  DPOWERX, 
and  DPOWERY  are  used  to  calculate 


f(x,y)  =  (x  +  iy)n 


Sx 


Sy 


respectively,  and  where  n  is  an  odd  positive 
integer.  The  entry  points  EPCWER,  DEPOWERX, 
and  DEPCWERY  are  used  when  n  is  an  even 
number . 


Function  subroutine  to  compute  (-l)n  where 
n  is  an  integer. 

Function  subroutine  to  compute  xn  where  n 
is  an  integer.  If  x  =  o  .  PWRF  will  return 
with  xn  =  o  ,  even  when  n  <  o  . 

Generates  the  normal  equations  by 
1  Transposing  matrix  A 

2.  Multiplying  and  forming  the  matrix 
product  A^' 

3.  Multiplying  and  forming  the  matrix 
product  A^b  for  the  right-hand 
side  of  equations 

Additional  Capability:  List  the  matrix  A  using 
the  subroutine  LISTARAY. 

A  subroutine  to  list  matrices,  labelling  rows 
and  columns . 


UCSD  Library  Subroutine.  Inverts  matrices, 
solves  linear  systems  of  equations,  and  com¬ 
putes  determinant.  It  uses  a  Gaussian  method 
with  two-dimensional  pivot  search. 
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ZEROLIST  UCSD  Library  routjne  will  zero  out  variables 
specified  by  the  list  defined  in  the  NAMELIST 
statement. 

PRTCOEF  Output  subroutine 

W1  Additional  entry  points:  W2  ,  TAUXZ1  ,  TAUXZ2  , 

TAUYZ1  ,  and  TAUYZ2  .  This  routine  computes 

wl  >  wI1  *  ^xz1  >  TxzIX  *  etc- 

GRID  This  routine  computes  coordinates  of  the  127  points 

used  in  the  finite  difference  method. 

SHEAR  Integrates  (Gaussian  integration)  t  along  right- 

hand  boundary. 

USAGE  OF  THE  PROGRAM 


Card  1 


Field  1:  VF  (Volumetric  Factor,  e.g  ,  0.6)  if  IFLAG  /  0 

or 

RHO  (Radius  of  Fiber)  if  IFLAG  -  0 

If  Vp  is  specified,  RHO  is  computed  as 


If  RHO  is  specified, 


Vp  is  computed  as 


Field  2 
Field  3 
Field  3 


E*  ,  Young's  modulus,  fiber,  psi 
V1  ,  Poisson's  ratio  fiber 
E**  ,  Young's  modulus,  psi,  matrix 
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.kpvm# 


Field  4:  v**  ,  Poisson's  ratio,  matrix 

Field  5:  IFLAG  ,  See  description  for  Field  1 

Card  2  Format  (14) 

IPRT  If  IPRT  *0,  do  not  print  the  matrix  of  the 

coefficients  for  the  linear  equations 

If  IPRT  <0,  print  matrix  of  the  coefficients 


Card  3 

Format  (14) 

NOMORE 

If  NOMORE 

<  0 

Call  GRID 

NOMORE 

«=  0 

Input  now  data 

NOMORE 

>  0 

STOP 

NOMORE  <0  calls  GRID,  which  causes  the  program  to  generate  the  coordi¬ 
nates  as  used  in  the  finite  difference  method  and  to  calculate  displacements 
and  stresses  at  these  points.  If  card  3  has  a  NOMORE  <0,  a  second  card, 

3a,  has  to  follow  which  has  a  NOMORE  2:0. 
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o  o  o  n'n  o  n  o  o  rVo  n|o  o  o  ono 


i  »  .  — - - - --- 

PROGRAM  ALBERTO  . .  . . . 

SHEAR  PROBLEM,  INTEGRATION  CONSTANTS  ARE  FOUND  BY  POSTULATING 
BOUNDARY" COWmOWS  AT'A^TTmTER- DF- TOTnYS ,  'SO  fHTAT'fOTAL- NUMBER  OF 
/  EQUATIONS  IS  MUCH  LARGER  THAN  NUMBER  OF  UNKNOWN 

COE  r  F I C I  ENT'S*."  EQU'ATTDnS*  ARE'  SOL  V'E'6"6y '  l'Ea'S  t"'SQUARE  S" METHOD 
NORMEO  GENERATES  THE  NORMAL  EQUATIONS _ _ 

***  IT  «  TOTAL  NUMBER  OF  EQUATIONS  GIVEN  BY  BOUNDARY  CONDITIONS 
“***"MBReO  V'nOmb'EK  Of  NORMAL  -  EQUATIONS  ( "NUMBER  0~F  UNKNOWN  COEFFICIENTS) 

'***  PROGRAMMED' SEPT  EMBER 'V3', *1967  . 

REVISED  SEPTEMBER  22,  1967 

.  MODIF  IED  0CT0Bfc'R'*9,  t96r  T0'uSE  *FUN'Crr£l'N';''rx  V  "lY)  A*(2*N+l ) 

_ MOO  I F 1  ED  OCTOBER  24,  1967  TO  USE  FUNCTION.  (X  ♦  IV)**N  IN  FIBE R_ 


M*  NOOD  •  MAXIMUM  ODD  POWER  USED  IN  BOTH  FIBER  AND  RESIN 

NEVEN*V*'M'A'jaM0M"fvE'N*lT0w'ER''USED*TN**Fi*6¥'ft  * 

- ‘CO  M  M  0  N  /  IN  PT/‘A‘ ,  BV  RH  OV  E 1TFWIY  *f  JVFNli  Y  . 

COMMON/MTKX/AAI  100,61),  RHS(IOO) 

T**w* - roCTT5  Th1T~ mATTmum . Nu'Hb ETOTT)01JNDlOfVTW)rfTWr"F'£'SMTT ‘f El) 

COMMGN/LABL/  INDEX21 100) ,  INDEX! 100) 

NAME  L  1ST  /  L  ABEL/I  N0E5f2  ,'TNDEx  . 

DIMENSION  BACK ( 100  I 

c  0  M  M  0  N/ N  ORM T'KX  /  AND  R'Mf  l'V6l' )"7'B'N  0  R  m'  ( 611  . 

COMMON/COE FF/A1 ( 20) ,  A2I20),  31120),  E01 

common/" nu'mbe r/nood ,  ufvtN,  nbr'eq,  "nva'r,  if 

COMMON/ PO I  NTS  'XX (150).YY(150),JFLAGI 150) 

0 1  ME  NS  I  ON’Tf'L  AGGT 5')T'J  F*0  R  M(*2  ) 

TYPE  INTEGER  STAR  _  __ 

.  DATA  f  IFL  "aG6~'«  ~8HF  1  B'ER'  ,  BHRESTn**  ',*  *  BHfNTRF  ACE , 

_ *  _ _ OH _ .  BH _ > _ __ 

UATA~(  JFORM~(  1  )  «  8H(1H»,95X) 

OATAINVAR  «  20),  ( NBREQ_«  61)t_LNC0L  LNUM_«  201 j _ 

" . *  . ( NO  D  D*"  1 9  )",*  ( NE  V E  N  ■  *20),  (IT  *  91) 

OAT  A ( PI  «  3.1415926)  _ 

6*  L  '■  LAPS TTmeTCLT 

_ CALL  START1HE _ 

READ  1,  A,E1,FNU1,E2,FNU2, IFLAG 

1  FORMATtSElO.4,110)  _  _  . 

. I  F  I  EOF,  50)  3lV  I",  4 

4  JF0RMI2J  •  8H_,A8t//)  _  . . 

C"" .  "" 

_ IF  I  1  FLAG)  11,  12 _ _ 

11  VF  ■  A 

RHO  ■  2  •  *SQRTF  ( 2.  *SQRTF(  3.  )*VF/PI ) _ _ 

. . GO  TO  13** . 

12  RHO  «  A  _  _ 

VF  *  RHO •  RHO ♦  PI /"( 8  •  * S QR f F*(  3*. l*f 

_ 1_3_GL  ■  El/12. ♦(  1.+FNU1)  ) _ 

G2 '  ■  E2/~{  2.  *  1 1  .+FNU2 ) ) 

G10VEKG2  «  G1/G7 
B  «  1.0 
A«B*SQRTF  I  j.  ) 

PRINT  2  ,'B",  VF', RH0,"E'iTF*NU1  ,*E2 ,  FNU2 Gl", G2 , G lOVE RG2 

_ L  FORMAT! 1H1, POX,  16H1NPUT  PARAMETERS ./ //20X , 1 6 ( 1 H- )/// _ 

18H  B  ,E20.5/ 

18H  VF  lE20.5/  . . . . 

28H  RHO  ,  E20. 5/ 

_ M.£l _ . . . 


noon  oor»o  r»oo 


c 


c 


48H 

NU1 

, E  20 . 5/ 

5BH 

£2 

#  E  20 . 5/ 

6BH 

NU2 

,£20.5/ 

IrtH 

Gl 

,£20. 

5/ 

1  8H 

G2 

,  E20. 

5/ 

18H 

G1/G2 

,£20. 

5//I 

LI 

*  L  / 

1000 

PRINT  603 ,  LI,  L2 
603  FORMAT  (SHOT  I  MT:,  I  13,  1H.  131 


L2  •  L  .  -  LI  *  1000 


PRINT  6C2  .  It,  NODO,  NCVEN,  NBREQ 

602  FORKATKIHmNUMRER  OF  BOUNDARY  EQUATIONS  SPECIFIED  * ,  I  5 , /26H0MA  X I  ML' 
*M  000  POKER  USED  «  ,  I  5 ,/ 26H0MAX  I  MUM  EVEN  POWER  USEO  » ,  I  5 , /20H0NUMR 
♦  ER  OF  COLUMNS  «  ,  IS) 


GENERATE  MATRIX  AND  RIGHT  HAND  SIDES  OF  EOUATIONS 


CALL  GETMTRX 


FORM  NORMAL  EOUATIONS 
CALL  NORMEO 


SOLVE  EQUATIONS  *  ' 

*"»  CALL  NWMAT INV( ANORM,  BNORM,  BACK,  61,“‘N0R£0Y  l.  DET,  IDE  T) 

S* 

I0ET20  *  20* I OEf 

PRINT  971.0ET, IOET20  _  _  _ 

971  format'  ( //V,  i'x7"Vh6eT*-,c20.sV5'x,  WhiTmTs  10**YI4) 
IFIOET  .NE.  0. I  GO  TO  318 

C 

317  READ  972,  N1 ,  STAR 

IF  (STAR  '.Eg.'  IH*)'6,"31'7  . . . . . 

972  F0RMATJ_A8j|_Al  1 _ 

c 

c 

C  «**  BACK  SURSTITUE  SOLUTIONS  IN  TO  "ORIGIN  ALMAT RTx . 

318  00  319  1*1  ,  I  T  .  .  . . 

BACK ( I  1  *  0. 

319  CONTINUE 


DO  320  JJ*1,  I  T 
DO  320  KK=  l , NBREO 

BACK  ( J  J  I  =  BACK  (  Jj'l  BNORMl KK )  *  A  A  (  jj ,KK ) 
320  CONTINUE 


:  ***  CALCULATE  SUM  OF  SQUARES  OF  RESIDUES . . 

SUMSURES  *  0. 

00  335  J  *  1  ,  I  T  . . . .  . 

_ 3 35  SUMSQRES  »  SUMSQRES  ♦  (RHS(J)  -  BACK(J)|**2 _ _ 

340  PRINT  904,SUMSQRES_  _ 

904'  FORMAT!  1  Hi",  BX.ThJ ,  UX'.ThBACK  ( J  jYlVx", 6HRHSI  J  )V25X',  1 7HSUM  RESIOUES  S 
*Q  ■ , E20. 5  •  /  )  _ 
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noon  n  o  'n  r»  Jo  o  no 


mmm mmmwvsm** 


PM  NT  90  3 1  (J,  (SACK  <J)(  i<HS(Jli  INDEXIJ),  J*1,IT) 
FORMAT  (110,  2E20.5,  IPX,  A8) _ 

STORE  CULFF I C I  ENT  S 


00  2  7  J* l , 00 


00  2R  JM.NVAR  _ 

A1  I  J  )  *  BN0R  K  t  j ) 

A2IJ)  -  BN0KM ( J ♦ NV  AR ) 
CONTINUE 

00  30  J«l,NUM 
B  l  ( J )  «  QNilRH  I  J  +  NCOL  l" 

CONTINUE  _ 

EOl  BNORM!  NBREU*) 


PRINT  THE  COEFFICIENTS 


CALL  PRTCOEF" 

call  shSar 


HEAD  5,  II 
IF  I  1 1 )  6,  AO,  AO 
5  FORMA  TWO* . 


AO  I  -  0 


99  I  •  !♦! 

X  >  XX(I) 

IFLAC  «  JFLAGi  n" 


V  •  Tflll 
. . LABELS 


INDEX?! I ) 


GO  TO  1 100 »  200,  100,  250 1  1  FLAG 

I  FLAG*  1  •V,Vl  BER*  . 

IFLAG-2... RES|N 
IFLAG*3... INI  ERF ACE* 

_ IF.  ,L  ApjfA  . . 

100  CONTINUE  _  _ 

PRINT  3,  I  FLAGG  1 1  FLAG),  X,  V 
3  FORMAT (/,1H0, AB,20X,  3HX  ■iE15.Al5Xi 3HY  »,E15.4> 

PRINT  JFORM,  LABELS . . . . 

_ CALL  WUX.Y.WWI) _ 

CALL  TAUXZl (X,Y»TAUXZZ1) 

CALL  TAUYZl ( X, Yj T AUXYV 1 1  _ _ 

PRINT  902 1 WMl , f AUXZ  Z 1 » TAUXVy  1 

902  FORMAT ( AH  HI* ,  El  5.  7,  SX , 10HTAU  XZ  l  *jE l 5. 7t5X, 10HTAU  VZ  1 

1  E15.7//J  . . . 

„G0  TO  199,  99,  201)  IFLAG . . 

200  CONTINUE  _ 

PRINT  3,  IFLAGGI IFLAG) ,  xV  ¥ . 

PRINT  JFORM,  LABELS  _ 

201  CALL  W2(X,Y,KW2) 

_  CALL  TAUXZ2IX.Y.TAUXZZ2) _ 

CALL  TAUYZ2! X, Y»TAUX YY2 ) 

PRINT  901,HH2,TAyXZZJ,TAUXYY2 

901  FORMA T ( AH  H2« ,E~15. 7, 5X,  10HTA0  XZ  2  *,E  15. 7, 5X,  10HT AU  YZ  2 

..-_.klhi.mi _ ; _ _ _ 
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o  n  n  n  n 


i 


I F  (  IFLAC  .Nfc.  3)  GO  TO  99 
C 

COSTllCTA  «  (  A-X)/i<HO .  V  S'lNTHETA  *  (IT-Yl/RHO 

FTAUt  *  TAUXZZl  *  COSTHLTA  ♦  TAUXYYt  *  MNTHF.  TA 
FTAu?  =  T AUX2/2  *  COSTHETA  TAUXYY2  *  S  INTMETA 

PH  I  NT  90b,  FTAUl,  FT  AU2 

90b  F  OHM  A  T  (24X,  10mF(  TAU  1)  615.  7,~5X ,  "l  OHFI  TAU  2)  -,E15.7,//I 

GO  TO  99  _  _ 

C 

250  HE  AO  5,  IFLAC  _ _ 

IF  (IFLAGI  TOO,  305,  315 

.  this'  VARIABLE  is  calleo  nomore  in 

IFLAC.  »  -1  CALL  _GRIO_ _ PI0.6RAM  DESCRIPTION 

IFLAC,  «  O'.  .  INPUT  NEW  DATA 

i FLAG  »♦!...  STOP 


300  CALL  ZEKOL I S T (LABEL )  t  CALL  GRID 

JFORM  (  2  )  -  8HYHY//T  »'  'GO  TO  'AO" . 

305  PRI NT906 

310  PRINT  500  »  GOTO  6 

906  FORMAT  ( ///,  1'HO  TiSi  ,‘JfOlT  C  ‘WO"  "  fOTN'T'SI 

500  FORMAT! 1HI) 


315  PRINT  906 

316  L  -"L'APsrrME  ai 


3161  LI  ■  L  /  1000  S  L2  -  L  -  LI  •  1000 

PRINT  AOSY'L'lT'L?  . . 


ENO 


j 

A 


1 


53 


SUBROUTINE  getmtrx 

C  ***  IN  THE  F  I  BEK  *  USE  THE  FUNCT_tON_ 

c  . “  IX  ♦  i v ) **m“  "V  h*“»  0717273', 

IN  THE  RFSIN,  USt  THE  FUNCTION 
C  (X  ♦  Iy;  **N  N'-'  r.3#5,' .  . 

C 

c 

COMMON/ INPT/AA,B,RHO»61,FNUl,F2»FNU2 

COMMON/ MTKX/  A ( 100,61)  ,'  KHSI 100) . 

NAMELIST/MATRIX/A,  RMS 

COMMON/NUMBER/NOOO,  NEVEN,  NBREQ,"  NVARi  I T 
COMMON/COMPLX/Z  < 20) 

tgui  valence  i  z.zevcN) .  . .  . " 

01  MENS  I  ON  ZEVEN ! 20 ) 

DIMENSION  XXI 20  ),YY  (20)VCT('?0)  ,ST(  20  f . 

COMMUN/POINTS/XPT (  ISO)  ,  YPTI1301,  JFLAGI 150) 
COMMON/LABL/LABEL  UbO')T  LA8EL2l  i'00') 

01  MENS  I  ON  SCR ( 100 ) 

01  MENS  I  ON'  NuMfiER  ( Ta")  . "*' 

OAT  A ( I H  »  62000000B),  (10  •  64000000B),  (IF  •  660000000), 

V  TIC  *  bTOOooooB) ,— m  »  Tioooo'oo'B'iTTl J  »  aioooooobi 

DATA! NUMBER  *  LIU,  1R2,  IR3,  IR4,  IRS,  1R6,  1R7,  1R8,  1R9 

"  *  . . *2ftl0 , 2RT77 2R  12, 2R  1 3",  2Rl4) . 

UATA(PI«3. 1415926531 

CALL  ZEROL I S  T (MATRIX)  t  CALL  INITIAL 

*  G  l «  E  l /  ( 2.V(T. ♦  fNGTf  I  t"G2»e2'/(2“.*Tl.  +  FNU2)V  *  CG*G2/Gl 
NCOL  •  2*NVAR 

NUM  ’•"NE  VEN . .  . .  . 

it 

C  EQUATION' V  (POINT* AT"W1-W2*0  . * 

c _ 


'  X« AA-RHQ "  i  Y-B 

c. .  .. 


2(11  •  XPT(l)  -  X  S 

LABEL ( 1 )  »  AKA 

2(2)  "•  YP T(l)  •  Y 

LABEL2UJ  «  BHA  1,2  -  W 

JFLAGI  1)  •  3 

11-2 

CALL  POWER 

c 

io 

00  10  N«l,NVAR 

A  ( 1 ,  N 1  •  2  ( N  1  S 

CONTINUE 

All, N+NVAR )  »  -2 ( N) 

ZEVENIll  -  X  $ 

CALL  EPOWER 

ZEVEN! 2 )  »  Y 

00  11  N-l.NUM 

A(l,NCOL*N)  •  ZEVEN ( N) 

11 

CONTINUE 

A ( 1  , NBREO )  ■  1. 

c 

c 

EQUATION  2  (POINT  A)  TAUXZ1/G1- 

TAUXZ2/G1-0 

c 

LABEL2I2 I  •  BHA 1 2  TAUX 

2(1)  ■  X  *  2(2)  ■  Y 

CALL  DPOWERX 

00  20  N«1,NVAR 

A ( 2 , N I  *  2 (N )  i 

A( 2, N*NVAR  )  ■  -2 (N)*CG 

20  CONTINUE 


5A 


oooooc  oor>  r>  non 


ZEVENI  11  *  X 
CALL  DEPflRcKX 
OU  21  N»  1  ,  NUM 
A ( 2  »  N  +  NCOL  )  »  ZEVEN(N) 
21  CONTINUE 


ZEVENI2)  «  Y 


EQUATION  3  (POINT  A)  TAUYZl/Gl-0 


LA0EL2 ( 3 )  «  8HAI  TAUY 
(III  *  X  %  2(2) 

CALL  OPOWERY 
00  30  N-l.NVAR  • 

A  (  3  i  N )  *  Z  ( N ) 

30  CONTINUE 


Z6VENI  1  >  =  X  : 

CALL  DEPOWERY 
00  31  N* 1  NUM 
A ( 3 , N+NCOL I  *  ZEVENIN) 
31  CONTINUE 


ZEVENI2I 


EQUATION  4  TAUYZ2/G2-0  IN  POINT  A 


LABEL2I4)  *  0HA2  TAUY 

2(1)  »  X  '  - 

CALL  OPOWERY 
DO  40  N« 1 iNVAH" _ 
A;4,N*NVAK)  «  Z ( N ) 

40  CONTINUE 


Y  . . 


EQUATIONS  5  THRU  18  POINTS  6  WI-W2-0 


GENERATE  THE  14  POINTS  XX  AND  YY 


DELTA-6. *PI/1B0. 


LABEL2 I  5 )  -  8H81.2  -  W _ 

DO  6  0  ~I  -  l « 1 4  i  FI*  I 

THE  T  A-F I *OELT  A  _ 

xx  ( i  )•  a  A-RHoicoVFVf  het‘a  i’T'c  t  (Ti -'cosfTthe  vaT 

VY  (  I  )  -U-RHO+S INF  (  THET  A )  *_STU  )  »S  INF  I  THETA ) 
CONTINUE  .  . . 


00  70  1-1,14 

Z(l)  »  XPTII  1)  -  XXII) 

LABEL  1  I  1  I  *''i  B '♦"nJmb'ERTI  )' 
JFLAG III)  »  3 

I1-U  +  i . 

CALL  POWER _ 

DO  60  N-l.NVAR 

A  I  l  +  4  ,*N )  *  Z  TnT . T 

CONTINUE 


Z  ( 2 )  »  YPTIIl)  »  YYII1 


At  I+4.N+NVAR)  ■  -ZIN) 


ZEVEJ1LI  I  -  XXII) _ i 

"CALL  EPOWER 
DO  61  N« l j  NUM 

a  ( i  n+ncolT'V  zevenTnT 
CONTINUE  _ 


ZEVENI2)  -  YY  1 1  ) 


A ( I *4 , NBR£  Q )  *  1. 

70  CONTINUt  _  _  _ 

C 

C  EQUATIONS  19  THRU  J2  f ( T AU1  )-F I T AU2 ) *0  POINTS  0 

C 

LABEL?  119)  -  8HB  £ (  T AU )'  _  _ 

00  90  l»l  ,14  . . . . . 

211  (  »  XX (11 _ % _ 2(2)  «  YY 1  I  1 _ 

. CALL  OPOWERX ' 

DO  71  N* 1 #  NVAR 
75  SCR  ( N)  «  7  ( N )  *  CTI  I  ) . 

c 

2EVENI  l  )  -  Xxfn" . ~i~  '  ZE YEN (  2  )  »  YVt  I  > 

CALL  DtPOwEnx _ 

“on  76  N»  1 .  NUM 

SCR(N*NVAK)  -  ZfcVEN(N)  *  CT(l) 

'  76  CONTINUE  . . . 

2(1)  *  XXII)  t  2(2)  ■  YY ( I ) 

CALL  OPOWERY . .  ' 

00  80  N-  1  , NVAR _  _ 

A  ( I  ♦  1 0 ,  N )“  «“SCR  fN )  ♦'  z7nT*S TTI  >' 

AlI^lS, N+NV AR )  «  -A ( I  ♦  1 8  ,  N )  ♦  CG 
80  CONTINUE  .  . 

c 

ZEVEN(l)  *  XX'O'V . "V  . . 2EVENI2  )""*  YY  (  I  ) 

CALL  OEPOWERY _ _ _ 

00  81  N* 1 , NUM 

A(  ( ♦1B,N*NC0L )  -  SCRIN+NVARl  ZEVENIN 1  ♦  ST « I ) 

■"81  CONTINUE" . . . .  .  "  . 

90  CONTINUE 

c  "" 

C _ EQUATI  ON  33  POINT  C  W1«K _ 

C 

X«AA  1  Y«B-RHQ _ _ _ 

C""~  " 

LA8EL2 ( 33 )  »  8HCi  -  W 

Z(ll  ■  XPTdi)"  «  X  i  Z(?)  «  YP Till) 

_ LABEL! 11)  -  4RC _ 

JF'.AGI  1 1 )  «  3 

11  ■  1 1  ♦  1  _  _ _  _ 

CALL  OWER 

00  I'.l  N»l,NVAR  _  _ 

A(33,N)  -  Z ( N ) 

100_C0NT  I  N'JE _ 

C 

2EVEN11)  «  X  »  ZEVENI2)  »  Y  _ 

CALL  EPOWER 

00  101  N-l.NUH  __  _ 

A(33,N*NCOL)  «  ZELENIN) . " 

_10I  CONTINUE _ _ 

"  “"A(33,NBREQ»  «  RHS ( 33 )  «  1. 

C 

c 

C  EQUATION  34  POINT  C  W2-K 

C'"' 

_ UmtiAAL  ?-.§HC2  ,  JZ—M. _ 

2(1)  ■  X  S  2(2)  «  Y 

CALL  POWER  _  _  _  _ 

DO  120  N» 1 ,  NVAR 

. AJMjJ&m&LjtJUJIi _ _ 


56 


UU  UOU  O  O  u  oo 


120  CONTINUE 

KHS ( 341  *  1. 


EQUATION  35  POINT  C  I  T AUYZ 1- T AUYZ2 ) 2G1 -0 
C 

LABEL?  (  35  )  »  tlHCl?  TAUY 

2(1)  *  X  *  2(2)  *  Y . 

call  upowfry  _ _ 

00  140  N  = 1 ,  NVAR 

A  (35,?;)  *  2  (  N )  S  A(  35,N*NVAR  )  •  -2(N)  *  CG 

140  CONTINUE 
C 

2EVENII)  *  X  S  .  2EVEN(2)"«Y  " 

CALL  DEPOn'ERY  _  _ _ 

00  141  N*l,NUK~ 

A(35,N*NCOL)  *  2E YEN ( N ) 

141  CONTINUE 

POINT  0  EQUATIONS  36  ThHU'44  TAU2Y1/GI-0 

GENERATE  THE  ~9’  COORD  IN  ATEIf  XX"  AND~YY  * 

XO  *  A4  -  KHO  i  "Y  *  B" 

1FIX0I  150,  150,  151 

150  I  DPT  *  12  " .  * - OPT  «~  lj; .  "»  GO  TO  152 

151  I  OPT  «  0  t  DPT  -  10. 

102  OELTAX  .RHO/OPT  "  ~  “ 

LABEL2I36)  ■  BMC'  TAUY" .  .  ‘ 

00  153  1=1, (OPT  l  FI  *  I 

XX  l  I  I  «XO*H  *OEL  T  AX  . . .  . . "" 

153  CONTINUE  _ 


00  160  1*1,1  OPT 

2(1)  *  XPT  ( I  I J  *  '  XX  I  I  ) .  T'  "  lUr-  Y PT  (  I  l  I  i-  Y 

LABEL! 1 1 )  *  10  ♦  NUMBER (  I) 

JFLAC(II)  *  1  . " . .  . 

II  *  II  ♦  1  _  _ _ _ _ 

CALL  OPOWERY 
DO  155  N=1,NVAK 

At  1  +  35, Nl  •  2  ( N )  .  . .  . 

155  CONTINUE  _ _ _ _ _ _ _ 

2EVENI 1 1  *  XX  (JJ _ S _ 2EVEN12I  -  Y _ 

CALL  DEPUWERY 
00  156  N* 1 , NUM 

A 1 1  ♦  35 ,  N+NCOL  f '•■■2W6N  ("nT .  *  ’  . 

156  CONTINUE 

160  CONTINUE .  “  ~ .  . 

_ 12  ■  36  ♦  IDPT _ 

EQUATION  45  POINT  E  TAUZY1/G1-0 _ _ _ _ 

X*  A  A  i  Y«B _ _ _ _ _ 

_ Z  ( 1 )  *  XPT  Ul>  ■  X _ * _ 2(2)  ■  YPTI  I  1 1  »  Y 

LABELTU)  «  4RE 

LABEL?! 12)  ■  BHE  -  TAUY _ 

JFLAG(llT*''l 
II  «  II  ♦  1 


57 


CALL  OPOWERY 
00  i/O  N«1,NVAR 
MI/.NI  -  Z !  N ) 
l  TO  CONTINUE 


ZEVFN(i)  *  X  $  ZEVENI2)  «  Y 

CALL  OE POWER Y  . 

00  171  N-1,NUM_  _  _ 

"  A  ( I  2  ,  n+ncolV*  ZfcVEN'lN  ) 

171  CONTINUE  _  _ 

12  «  12  ♦  I . 

c  _ 

c  EUUATION  46  POINT  E  Wl-K 


LABEL2II2)  *  8HE  -  Wl 
Zll)  *  X  i  2(2)  «  Y 

CALL  POWER . 

00  190  N»1 ,NVAR 
A  (  1 2  ,  N )  i  'i'fNf 
190  CONTINUE 


ZEVEN(l)  -  X  S 

’  CALL  EPOWEff . 

OU  191  N»1,NUM 
A 1 1 2 ,  N+ NC'OLT*  «"7"EVE  ti  ( NT 
191  CONTINUE 

A  t  1 2 ,  N BR'fi)  I  “  R h'sTT 2)  *  1 7  " 


ZEVENI2)  « 


EQUATION  S"47'  "Th  R~U  '  55"  PO" IN' f  5  F  wW 


L  A8  E  L2  Cl  2  + 1 ')”'»”  8H  F 

X-AA  i  YO-R  t  DELTA 


IFPT  -11-1 
00  200  1*1,9  *  FI«i 
YYm-YO  +  FI*DE'LTAY 


8HF  -  W 
OELTAY— RHO/IO 


DO  220  1-1,9 
ill)  »  XPT(Il)  - 
"LABEL!  ID  ■  IF  ♦ 
JFLAG (II)  •  l 

u  -  u  ♦  i " ' 

CALL  POWER 

oo  210  n»T,War" 

A ( !♦ 1 2,N )  -  ZIN) 


CONTINUE 


YPT(ll)  -  YY I  I ) 


NUMBER! I ) 


ZEVEN(l)  -  X  » 

CALL  EPOWER 

00  211  N«l,NUM . . . 

A  ( I  ♦ 1 2 , N+NCOL )  -  ZEVEN(N) 


CONTINUE 

Aim?,  NBRfcQ  I  -  RHS  (  l  ♦  I  2  )  ■  l. 

CONTINUE  . . 

12  -  12  *  9 


ZEVEN! 2) 


XPT 111)  -  AA 
LABEL!  ID  •  ARC 
JFLAG (II)  •  3 
It  -  II  ♦  1 


YPT(Il)  -  B-RHO 
ICPT  -  I  l  -  l 


000  o  O  O  J  i  o  o  o  o 


C  EQUATIONS  56  THRU  59  POINTS  G  W2*> 
C 

LAI3EL2 (12*1)  »*'8HG  -  '  w'2 

YO-U-RHO  i  OELTAY— t?.*B-KH0l/5. 


DO  250  1*1(4  i  F  I  *  I 
Y«YO  +  FI*f)ELTAY 

JFLAGUl)  *  2  _ 

LABfcL(U)  *  IG  ♦~NUMflER(  I  ) 

c 

ZU)  *  XPT  U  l )  *  X .  t . . Z  (21  "»  YPT(ll)  ■«  Y 

11  *  I  1  ♦  1 
CALL  POWER 

_00  240  N* 1 ,NVAR  _ 

AT  I  *  1 2 iN  +  NVAR )  *“ Z~( N ) 

240  CONTINUE 

RHS  (  I  *  1 2 )  *  1.  '  .  . . . 

250  CONTINUE 

12  *  12  ♦  5 

H 

EGUA  T I  ON  60—  POI  NT- H-W2*K  ~  ' 

X«AA  i  Y-'-B  .  "  . 

LABEL (II)  »  4RH 

LABEL?  ( 12  )  . * . .  * . 

JFLAGUl)  *  2 


Z(l)  *  XPT(Il)  -  X  t  Z 1 2 )  «  YPT(ll)  *  Y 

11  *  II'V'T .  . ~ . 

CALL  POWER 

00  260  N*T#fiVAir' . 

A (  I  2  i  N+NVAR)  »  Z(N) 

Ybir CONtTnu? 

RHS  (12)  -  1.  _ _  _ 

12  -  i2  *v*i  .  . 

EQUATION  6 1  P  0 1  N T"*H  f  A  U  Y  2 2V 6*2 " 0  . 

L AOELZTt 2 )  *'  8HH  -  TAUY 
2(1)  •  X  i  2(2)  «  V 

CALL  OPOWfcRY' . . . . . . 

DO  280  N* I ( NV AR 

AT  1 2  » N+NVAR  )  »  Z  (  N) . . . . 

280_C0NT I NUfc _ 

I 

EQUATIONS  62  THRU  73  POINTS  I  TAUYZ2/G2-0 


XO  »  AA  -  RHO 

I F  (  XO )  ?90,  '"90  r  "291 . . 

290  AO  «  XO  ♦  AA _ t  XO  »  -XO  t  IPT  »  0  A  GO  TO  292 

291  AO  -  A A  t  XO  •  0.  »  IPT  •  1 

292  OELTAX  •  AO  /  13.  _»  __Y  *  ^.B _ 

. L  ABEL2  (  I  2*  iT  '*"  8H  f '-""f  ALIY  . 


00  300  1-1,  12  »  F I -  I 

X  «  XO  ♦  F I *OEL  TAX _ 

JFLAGUl)  -2 

LABELUl)  ■  I  I  ♦  NUMBER (  I  ) 

z  rrr-'xp'TY  n  rrx  r  zT2T"-'ypt'(  it  > v 

ii-ii*i 


59 


°  °°,  |  °  I  !  on  o.o  o,  o  non 


CALL  UPOWLKY 


1)0  245  N*  l , NVAR 

All* I2»N*NVAR)  «  Z ( N I  „  _  _ 

295  CONTINUE 
100  CUNT  I  HUE 

IF(IPT)  301.  321  .  .  . 

tOUAT  1  ONS  7 6  ~T  HR  U  76  PU f N"fS~ j~T AtjVz 2Tc2*0 

301  DEL  T  AX  «  (AA-RHUI/3.  V  Y  «"  R 

LA0EL2I76)  »  0HJ  -  TAUY 


DO  320  I  -  1  ,  3_  t  F  1»  1 _ 

x«  (  fi  -n  *otL  Tax  *  .  ' 

JFLAGII1I  «  2 

LABEL  (III  • '"fj'T  NUNOErViV . " 

2(11  *  XPTdll  «  X  %  212)  »  YPTIIl)  »  Y 

II  ■  I  1  ♦  1™ .  . . 

CALL  OPOWERY 


00  310  N* 1 jNVAH 

A  ( I  ♦  7  3  »  N*NVARf~  *  *2  ( N  I  . 

310  CONTINUE 

320  CONTINUE . 

"•♦•"EQUAT  ION- 77-B6~bW f/DYV0"ArtfN6'  POTNT STANDI  N  POINT  C 


321  LABEL2I77)  »  flHF  OW/OY 

*“  ‘00  160  I  « 1,10  . . 

_ 2(1)  •  XPT ( 1 F PT* 1  I  S  2(21  «  YPT(IFPT*I) _ 

CALL  OPOWERY 

00  330  N-l.NVAR  _ _ _ 

A  ( I  ♦  76 ,  N*l  •  2'(N) 

.330  CONTINUE _  ...  _ 

_ 2fcVEN(J  I  «  XPT  (  1FPT  *  I  I _ | _ ?EVE^;2)  »  YPT(IFPT*I) 

CALL  DE POWER Y 

00  331  N« 1  * NUN  _  _ 

A ( I* 76 . N+NCOL  I  *  2EVEN(N) 

131  CONTINUE  _  _ 

3*0  CONTINUE 

_ LABE_L2_(  86  I  ■  8HC1  OW/UY _ _ 

•  ••  EQUATIONS.  87  .- .91.  OW2/OY*O..IN  .POINT..C  AND  ALONG  POINTS  G 

L A8EL2 ( 87 )  -  8HC2  OW/OY 

X*  AA  $  YO-B-RHO  t  OELT A Y»- f'2 . •B-RHT )> 5. 

_D0  650  1  ■  1  .5  S  FI  «  I _ _ 

2(11  •  XPTlfCPTM)  *  2121  -  YPT<  fCPT*l ) 

CALL  DPOwERY  _ 

DO  660  N«1,NVAR' 

A ( I  *  86 i N*NV  AR )  ■  2(N) 

660  CONTINUE' 

650  C ON TJNUJ _ 

LABEL2I88)  *  8HG  OW/OY 

JFLAGJUl  ■  *  . . . . . . . .  -  ... 

RETURN  t  END 


60 


SUBROUTINE  PURER 


c 

C  **♦»  This  SUBROUTINE  COMPUTES  VALUES  OF  THE  FUNCTION 

(X  ♦  IY)**N  ,  N  »  1,  2,  3,  .  .  . 

C  ANO  ITS  FIRST  DERIVATIVES  WITH  RESPECT  TO  X  AND  Y 

C 

CUMMON/CUMPLX/Z ( 10) 

CUMMUN/NUMBCR/NUDO,  NEVCN,_NBREO,_  NVARj_J  T 
TYPE  COMPLEX  2,21 
EUUIVALENCF  (Z,  OZ  ) 

DIMENSION  UZIPO),  CI230) 

COMMUN/LABl/ AI2001 
NAMEL I ST/LABEL/A.C 
C 

IF  I  Z  I  1  )  .  EG.  fo.V’o.  n  GO  TO  no" 

Zl  *  2(11  *»2  1  •  NUM  •  (NUD0*l)/2 

GO  TO  10 
C 

C  A** 

ENTRY  EPOWER _ 

IFIZlil  .GO.  (0.,  0.1)  GO  10  no 

Zl  •  zm  **2  .  i  "  NUM  "*  NGVEN/2 

2(1)  «  Zl 

C  - - -  ” - - - - - - 

10  00  20  I  * 2 , NUM 

2(1)  *  Z7TT)~*"Yl 

20  CONTINUE  A  RETURN 

-  **• 

ENTRY  DPUWERX . . 

c 

X  *  02(1)  S  Y~  *  Ol  f?  ) 

□2(1)  »  1.  A  02(21  *  0. 

IFIX  . eo.  o. ""  Vano." . "y"‘.eo."  o.T  go"to  i in 


NO  *  0  A 

N 1  - 

3  A 

N  2  - 

NOUO 

**• 

GO  TO  50 

ENTRY  DEPOWEKX 

X  •  02(1) 

IFIX  .EO.  0. 

A 

•  AND. 

Y  .EO. 

Y  •  OH?) 

0.)  GO  TO  110 

NO  •  1  A 

Nl  * 

2  A 

NZ  - 

NEVEN 

50 

00  70  N«N1 ,N2,2 

NCOMB  *  N  *  (NMI/2  - 
DO  65  I *1 ,2  1 

2 

i r*  t-i 

1 

NN  »  N- 1 

SUM  »  0. 

55  DO  60  K«IT,NNV2 .  i" .  frV'N-K 


SUM«SUM*  FK  *  C ( NCOMB ♦ K  I  *  PWRF  ( X ,N-K- 1 )  •  PWRFIY.K)  •  AMINUSIK/21 
60  CONTINUE" 

C 

65  OZ  ( N«Tl -NO)  *  SUM .  . 

70  CONIINUE 

RETURN  .  .  .  . " .  . . 

C 

c  *•"• 

ENTRY  OPOWERY 
C 

X  -  02(1)  A  Y  »  02(2) 


61 


r»  n  o  o  r>  »  o  o 


om  >  ■  o.  »  uz(2i  *  *i. 

I F  {  X  .FO.  0.  .AMI).  Y  .60.  0.)  GO  TO  110 

NO  »  0  4  N  f  ■  3  "»  '  N2  «  NODQ  ‘ 

GO  TO  80 
C 

tNTKY  DtPOWERY 
C  ‘ 

X  *  02(1) _ i _  __  _ Y  »_02I2)  _ 

I F  (  X  .60.  0.  ._An6.  Y'  .60.  0V>  GCl'TO  110 

NO  *  1  V  N 1  *  2  4  N2  -  NtVEN 

C 

80  UO  100  N«N1 .N2i  2 

NCOMU  «  N  *  IN* l 1/2  -  2  ‘  1  NN  *  N 

_U0  90  1-1,2  4 _ II  «  1-1 _ 4. _ SUN  *  0. 

C 

85  DO  90  K- I  I  ,NN,2  4  FK  *  X 

SUM-SUM*  FK"*C( NCOMBVkV  V'PWRF ( X,N-K)  *  PWRF(V,K-1) 
90  CONTINUE 
C 

_ 95_D2 ( N*I 1-NO)  »  SUM _ 

lO'O  CONTINUE 
RETURN 

c"' 

U0  00  120  N-3,20 

02  (N)  *  0."  . 

120  CONTINUE _ 4 _ RETURN _ 


•  •• 

ENTRY  INITIAL 


_ _ CALL  7ER0L 1ST (LABEL  1 _ 

A  I  2  )  •  A  I  3~)  «  1  . 

NMAX  «  XMAXOF(NO_OOt  NEVEN1  ♦  3 
IPOS  »  1 

DO  150  N« 5. NMAX  $*  NN"-  N  -  2 

_ _ 00  145  J- 1 ,  NN _ 4 _ K  »  N  -  J _ 

AIK)  «  AIK)*  A ( K-  1  I 

I  F  I  A  (  K  1)  1 45  f„14_5,_  J44 _ _ 

THE  ARRAY  C .CONTAINS.  THE,  .BINOMIAL.  COEFF  IC  I  ENT  S 
l  A4 C  (  I  POS  I  AIKJ _ _____ 

ipos  «  ipos  ♦  i 


145  C0NTINUF 
150  CONTINUE 

PETURN  "$ . . . fcND 


AMI NUS ( K/2 ) 


62 


GO'  TO  15 


FUNCTION  AMINUS(N) 

I  F  (  N )  5,  25,  10_ 

5  M  «  -N 
10  H  «  N 

15M*  M  .AND.  I  »  . IF  (M)  20,  25 

20  AMINUS  ■  -1  i  RETURN 

25  AMINUS  *  1  %  RETURN'  '  %  END 


10 

20 

30 

40 

50 

60 

70 

80 

90 


FUNCTION  PWRE (X,N) 

PWRF  *  1. 

IF  (Nl  10,  90 

ifixi  to,  ?rr 

PwRF  *  0.  I  RETURN 

IF ( X-  1  • I  AO,  90 
IFIN)  50,  90,  70 

M  -  -N  i .  no  60  1*1  ,M 

PWRF  *  PWRF  /X  »  RETURN 

00  80  I  *  1 , N  ~  ~ 

PWRF  *  PWRF  *  X 

RETURN  . . .  . . .  t  NO 


:0  o  o  oor>  o  n  n  o  <■.*-»  oo 


SUBkOUUNE  nokmeo 

FORMS  THE  NORMAL  EUUAT IONS  _F  OR  THE  SOLUTION  OF  THE 
OVEROE  TEKMINKD  SYSTtM  “  '  .  . 

COMMON/NORM TUX/ ANORM  (61,61) ,  BNORMt  61  ) 

COMMON/ MTRX/AAI  100,61)  ,  RHSI 100 ) 

COMMON /NUMBcR/NODO,  N  E  V  6  N ,  ’ N BRE Q ,  NVAR,'  I  T 
DIMENSION  ATR(61,100J _ 

00  20  J  «  UNBRED 
DO  10  K  >  1 , NBREO 
10  ANORM ( J , K)  •  0.0 
20  BNORMIJ)  •  6.0 

""TRANSPO srfH'f REOUCTO  mTtRI  X  . 

ISO  00  200  J«l,Tf' 

00  200  K  •  L.NHREO 

ATR(K.J)  ■  AAIJ.KI . . .  . . 

200  CONTINUE 


MULTIPLY  ATR  TIMES  AA 
00  360  J  •  I ,  NBRfcO 

OU  340  K  «"Y,NBR£Q  .  . . 

UO  330  KK  «  1, I T 

330 "ANORM I  J « K)  ■  ANORM!  J_,K  I  ♦  AYRTj •  KK  )  *AAfKK , K) 
340  CONTINUE 
3 SO  CONTINUE 

FORM  RIGHT  HAND‘YrO"£f'0¥'MRMAC‘'feOU'AI!dNS 
00  380  j“«"  iTnBRTO 

00  380  K  ■  1 , 1 T  _  _ 

380  BNORMIJ)  .  BNORM  ( J  )"”♦  ATRIJ.tO  •  RHSTk') . 

110  READ  TOO  ,1  lV  "i  <?" . "" .  . 

TOO  F OR M  A  T I  2  1 4] _ 

(Fill)  390,  400,  400 

390  CALL  L  l  ST  ARAY  |  A  A 1 100 ,  J.T  i/i?  R.6  0 1  JL?_ . 

400  RETURN  i  END 


I 


SUrtRUUM  NE  NWMAT  JNV! A,B,  [ NOE X , NMAX , N . K , DE T ERM ,  lOtT) 

NWMV 

1 

D  I  Ml \S  I  ON  A ( NMAX,  1 ) , R<  NMAX, l ) , INDEX (  1 1 

NWMV 

2 

EUUI  VALENCE!  IKOW, JROWVlRl,  (  rCOLUM.'jCOCUM,  1C  1 

NWMV 

3 

EOUl VALENCE!  AM AX, T, SWAP,  1  AM A  X ) , (PIVOT, TEMP, l TEMP) 

NWMV 

4 

OATA  (M|NUS*6000000000000000SI 

NWMV 

5 

c** 

ini  riAii/ArioN. 

DE  T  E  RM* l . 0 

NWMV 

b 

I  U6  T  «  0 

00  20  J*  l  ,N 

NWMV 

7 

20 

INDEX ( J ) 'MINUS 

NwVv 

8 

00  560  I » 1 ,N 

NwMV 

9 

C** 

SEAKCH  FOR  FLFMENT-  OF  LARGEST  MAGNITUDE. 

AMAX'0.0 

NWMV 

TO 

00  106  J1*1,N 

NwMV 

l  1 

IF (- INDEX! Jl )  )  106,  105,60 

NWMV 

12 

60 

DO  100  K 1  *  1  ,  N 

NWMV 

13 

IF (-INDEX  IK  l  )  )  100,100,80 

NWMV 

14 

80 

TFMP=A(K1 , J1 ) 

NwMV 

15 

IFITFMPIh}, 100, 82 

NWMV 

lb 

82 

TEMP.-rFKP 

NWMV 

17 

83 

IK-1  TtMP-IAMAxl  100,  100,84 

NWMV 

18 

84 

AMAX*-fCMP 

NWMV 

19 

1R«K  l 

NWMV 

20 

I  C»  J 1 

NWMV 

21 

100 

CON  f I NUE  - 

NWMV 

22 

105 

CONI INUE 

NWMV 

23 

"  “ 

I  F  (  AMAX  >120,115  . .  " 

NwMV 

24 

115 

DEIEKM'O 

NWMV 

25 

I DE  r  «  0 

RETURN 

NWMV 

2b 

120 

1  R0W  *  I R 

NWMV 

27 

ICOLUM* IC 

NWMV 

28 

INDEX!  ICOLUM  I  *  I  NDEX  (1  COlUM  I  .  AND.  .  NOl  ."Ml  NU  S 

NWMV 

29 

I F ( • NOT • (  IROW-ICOLUM) 1260, 140 

NWMV 

30 

140 

UETFRM.-OFTERM 

NWMV 

31 

C** 

EXCHANGE  ROMS. 

00  200  1*1 ,N 

NWMV 

32 

S  W  A  P  ■  A  !  I KOU.L  ) 

NWMV 

33 

A(iRUW,LI*A(  ICOLUM, LI 

NWMV 

34 

200 

A  (  ICOLUM,  U*  SWAP 

NWMV 

35 

I F ( .NOT.M)2o0,210 

NWMV 

36 

210 

00  260  1*1,  M 

NwMV 

37 

SWAP*  R(  1  ROW, L 1 

NWMV 

38 

RllROW,l)*R( ICOLUM, LI 

NWMV 

39 

250 

B(ICOLUM,L)*SWAP 

»’VMV 

40 

C*» 

SAVE  PIVUT  INFORMATION.  00  DETERMINANT. 

260 

INDCXI  I  I  *  !  ROW*  1 000000008*  I  CGlTjM*  I  ND^Xl  I  ) 

NWMV 

41 

PIVOT  *A(  ICOLUM,  ICOI.UMI 

NWMV 

42 

0ET£KM»D6TERM*PIVUT 

NWMV 

43 

c 

C"  •* 

KEEP  DETERMINANT  BEfwEEN  1.0E-20  AND'  l.OE-20 

- 

270 

DETERM1  *  ABSF(DETEHM) 

I F ( DE  TE  RM  1  . f  T  .  1.E-20I  GO  TO  275 

OETtRM  -  uETERM  •  1.E20 

10ET  *  IUET  -  1 

GO  TO  270 

275 

I F ( OE  TERM  1  .LT.  1.E20)  GO  TO  300 

DETERM  *  OETERM  /  1.E20 

I OE  T  «  IUET  ♦  l 

GO  TO  270 

65 


1 


300  CONTINUE 

C**  REDUCE  LEADING  COEF  TO  1. 

■  ~a  i  icolum  nmuTo  -TTO 

DO  350  L-l.N 

350  A(  ICOLUM, L')*HYC0XuM7lT/R"|WT . 

IFI.N0T.M13B0.360 

'  360  00  370  L-l.M . 

3  70  B I ICOLUM, L  1 »  B  (  1 COLUM , L ) / P I VOT _ 

C**  "SUBSTITUTE  FOR  NTH  VARIABLE". 

380  DO  550  L  1*  l , N 

I  F  I . NOT . I L 1-T COLUM )T~5 50 ,400~ . 

400  T*A(L1 , ICOLUM) 

AIL  1,1 COLUM) *0.0  . 

_D0  4  SOL*  1  ,N _ _ 

450  A(Li,L)*ATLl,L)-AUTOLUM,L>*T 
IFI.NUT.M)550,460 

460  DU  500  'LViVM"  . 

500  B(Ll,L)*B(Ll,L) -A (ICOLUM, L)*T 

550  CONTINUE"' . . . . . . 

C**  UNOO  ROM  EXCHANGES. 

L*N 

DO  710  L2-1.N 

jrow* i noex (l ) /Todooooooft  . . 

JROW* JROW. AND. 77777B 

JCOIUM* I  NOEX  (  L  )VAND .  77  7  7  7  A” .  . . .  " 

IFI .NOT. I JROW- JCOLUMI >710,630 _ 

630  DO  705  K»  1~»N 
SWAP*  A ( K , JROW ) 

AIK , JROW  >• A I  K i J  C  0 L  U H 1 . 

705  UK, JCOLUM)*SWAP 

710L-L-1  . . . . . 

740  RETURN 
END 


NWMV  44 
NWMV  45 
NWMV  46 
NWMV  47 
NWMV  48 
NWMV  49 

NWMV  50 
NWMV  51 
NWMV  52 
NWMV  53 
NWMV  54 
NWMV  55 
NWMV  56 
NWMV  57 
NWMV  58 
NWMV  59 

N«MV  60 
NWMV  61 
NWMV  62 
NWMV  63 
NWMV  64 
NWMV  65 
NWMV  6 t 
NWMV  67 
NWMV  6e 
NWMV  65 
NWMV  7C 
NWMV  71 
NWMV  "72 


•T 


SUBROUTINE  LIST  AR  AY  (A,  NMAX,  M  ■  N*"  1ST  ART ) 

C  ___  _  __ 

c.  ***  SET  UTDEOU'  "WHcN  U SI N C  Y H f S  SUBROUTINE 


FOR  USE  With  PROGRAM  ALBERTO 

DIMENSION  AINMAX 
DIMENSION  1 C I  10) 
OATA  (  (FORM  • 

,  11 

,  I  FORK (10) 

BHdHll  , 

BHl/.brtfattO, 

(1H  , 

8HH  COL,  I, 

BH4 , 3X ) , faH, 

8H  ROW)  , 

SHI  14,  IX,  , 

8H  , 

8H!fa)  , 

2 

I  FORM ( 3 1  * 

8HI1H5I  1 

8HW  ,10(5 

IFORM(n)  • 

00  10  I  *  1  ST  AR  T  , 

8H10E12.3, 

N,  10 

IM1  .  1-1  l 

IF  I IP9  .GT.  Nl  fa 

IP9- ~ f*9  t  LAST  ■  10 

,7 

b 

900 

IP9  •  N  »  J  •  N- 1 S lAR  T ♦  l  i  LAST  *  J- < l 0* ( J X l 0 ) I 

ENCOOt ( 0 1  900,  I FORM (3)1  LAST 

ENCODE  ( 8,  901V  7F0RM I d ) l  'L AST  . . . . . 

FORMAT  ( 4HW  ,,12,2H<5) 

901 

7 

8 

FORMAT  I  12  , b HE  12. 3,')  ' 

DO  «  I CO'JNT  •  1,  LAST 

1 C ( (COUNT  I  «  I Ml  ♦  I COUNT 

PRINT  1  FORM ( 1 ) 

PRINT  I  FORM!  2 )  , 
PRINT  I  FORM!  7 1 • 

riClJl,  J-ULAST1 

(J,  ( A  I J , I I ) ,  I  I«I,  IP9),  J,  J»1,MI 

10 

PRINT  I  F7)Rm(  2  1  , 
CONTINUE 

T I CTT )  J«l,LASfl 

PRINT  I  FORM  (10) 
RETURN  i 

END 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

lARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 

LARY 


3 

4 

5 

6 
T 

a 

9 

10 
1 1 
12 

13 

14 

15 
lb 
17' 
IB 

19 

20 
21 
22 

23 

24 

25 
2b 
27 
28_ 

29 

30 

32 


f 
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SUORUUI I NE  PRTCOEF 

COHMON/COEF  F  /  4 1  ( 20 )  V  42  (>0')T  81  <  20')T  EO 1  ' 
COMMON/ NUMflER/NOOO,  NEVEN,  NBREQ#  NVARt  IT 

C  ” 

PRINT  901 

901  FORMA T  (I H 1 , 2 IX, 2H»l‘,'20X ,>H8lV30xV2H»Yr ////  f  " 

C  _ 

NKAX  »  XM AX 6V t  NO 00 ♦  !  V  ~ NE Vic  N)~ . 

c _ _ _ _ 

PRINT  902  ~IN,  A 1 ( N ) ,  81IN),  A2IN),  N«1,NMAXI 

902  FORMATmOi2E20.<i,  10X,  E20.4I  _ 


PRINT  90  J •  E01 

901F0HMAT  (  i  f)  .'i>H  ‘For*  .  E2hV<. .  /Vi KiT 


SUBROUTINE  Wl(X,  Y,'  ANSI 

COMMON/  1NPT/A,B.RH0.  El ,  FNU 1 , E2 , FNU2 


C0MM0N/C0FFF/A1 120  )  ,  AYT20  f7~B  lTISl ,  £1)1 
COMMON/C CMPLX/Z (201 

c 

DIMENSION  ZEVENI20) 

EQUI VALENCEIZiZEVEN) 

COMMON/ NUMBER/ NOOO (  NEVEN,  NftREQ,  NVA'R,  IT 

. 

NZEKO  «  NOOO  ♦  1  %  NUM  ■ 
SUM  *  0 .  i  III)  >  X  A  Z 1 2 ) 
CALL  POWER 

00  10  N* 1 » NZ  ERO 

NEVEN 
»  Y 

• 

10 

SUM  »  SUM  ♦  A 1 ( N 1  *  ZlNi 

CONTINUE 

c 

ZEVEN(l)  *  X  A  Z  EVEN ( 2  1  ■  Y 

CALL  E POWER 

00  11  N-l.NUM 

li 

SUM  *  SUM  ♦  BUN)  ♦"ZEVfcNIN 1 ' 

CONTINUE 

-  c 

ANS  *  SUM  ♦  E01 

* 

c 

RETURN 

c 

ENTRY  W2 

C" 

NZERO  *  NOOO  ♦  l 

SUM  *  0.  5  Zfn  «  if  i  2(2) 

CALL  POWER 

vy 

DO  20  N* 1 i Nz£RO 

SUM  *  SUM  ♦  A2INI  ♦  ZIN1 

c 

20 

CONflNUT 

ANS  «  SUM 

RETURN 

c 

c 

c 

ENTRY  TAUXZ1 

NZERU  *  NODU  *  1  *  NUM  « 

GI*E1/(2.*(1.  +  FNU1!  1 

NEVEN 

SUM  .0.  $  Z  (  1  )  *  X  t  z  (i  1 

CALL  OPOWERX 

■  V 

00  30  N*l, NZERO 

SUM  ■  SUM  ♦  A 1 ( N 1  *  Z(N) 

c 

30 

CONTINUE 

/EVEN! 1 )  «  X  i  ZEVENI2)  •  Y 

Call  oepuwekx 

00  31  N* TTNUM 

SUM  «  SUM  ♦  Bl(N)  *  ZEVENINl 


*  3) "’CONT  I  NUE  *  . 

r 

ANS  ■  Gl'  V'SuM 

_ RETURN _ 

C 

c 

. ENTRY'"!  AU  XZ  2 

c _ 


Will  imnmnr^., 


C2»E 2/  ( 2.  ♦  l 1 7+FNU2 "if" 

NZERO  ■  NUOO  ♦  1 

SUM  >0.  »  2(1)  -  X 

CALL  OPOWERX 

t 

2(2) 

■  Y 

00  40  N*  1  *  NZERO 

SUM  -  SUM  ♦  A2(N)  *  Z  <  N ) 

C 

40  CONTINUE 

ANS  «  C 2  ♦  SUM 

RETURN 

C 

C 

*♦* 

C 

ENTRY  TAUYZl 

Cl  *  El / ( 2 .* ( 1 • +FNU1 ) ) 

NZERO  -  NOOD  ♦  1  S 

NUM  « 

NEVEN 

SUM  ■  0.  i  Z(l)  ■  X 

CALL  OPOWERY 

i 

Z  (2) 

«  Y 

DO  50  N»ltN2ER0 

SUM  »  SUM  +  A 1 ( N)  *  2(N) 

n 

«# 

50  CONTINUE 

ZEVEN(l)  >  X  J  Z EVEN ( 2 1 

CALL  OEPOWERY 

a 

Y 

00  51  N.'lTNuM'  . - .  . 

_ SUM  *  SUM  ♦  BUN)  ♦  ZEVEN(N) _ 

51  CONTINUE 

c  _ 

. ANS  •  *Gl  *  SUM' . .  ' 

RETURN 

"c" . . .  . 

C  *** 

ENTRY  TAUYZ2 

C  _  _  _ 

C2«E2/(2.*(i.*FNU2) ) 

N2ER0  -  NODO  _♦  1  _ _ 

SUM  *  0.  '*  2(11  -  X  *  2(2)  «  Y 

_ C A L L  OPOWERY _ 

00  60  N»1 ,N2ERO 

SUM  «  SUM  ♦  A2 ( N )  *  Z  <  NJ _ _ 

60  CONTINUE 

c  . . . 

ANS  »  G2  ♦  SUM 

c _ _ _ 

RETURN  t  ENO 


non 


SUBROUTINE  SMEAR 

COMMON/ I  NPT/  A,  R,  RHO,  El.  FNU1.  E2,  FNU2 
OIMtNSION  Ail'S  f.'  Xl(5'l,  A4I20),'  X4I20') 


DATA  (Xl« 

GLO 

8 

0 

.97)9065285, 

'.8650633667,  .  6794095683,  '.‘433  395394 1  ,  .  1480743390) 

GLO 

9 

DATA  ( A 1  = 

GLO 

10 

0 

•  0666  71  3443  , 

.  1494513491 ,  "  .'2 190  86  3625  ,  .2692667193,  .29552422471 

GLO 

12 

DATA  ( X4  = 

GLO 

19 

0 

.990237  7097, 

~  990 7  26  2  30  7  ,  .97  72  599 5007_.’95  7 9 168  192  ,  .9320  12  8083, 

GLO 

20 

u 

.9020906070, 

.8659595032,  .0246122308,  .7783056514,  .7273182532, 

GLO 

21 

0 

.671 9966846 , 

. 6125538897, "^5494671251 ,'  .4830758017,  . 4  1  3 77920‘. -  , 

GLO 

22 

0 

.3419940903, 

.2681521850,  .1926975807,  .1160840707,  .0387724175) 

GLO 

23 

DATA  ( A4  = 

GLO 

24 

0 

.00432  12771  , 

.0104982845,  .0164210584,  .0222458492,  . 02 793 700 ,0 , 

GLO 

25 

0 

•03)4601953, 

.0387821600,  .0436709082,  .0486958076,  .0532278470, 

GLO 

26 

0 

.03  74)97691  , 

.0613062425,  .0648040135,  .0679120458,  .0706116474, 

GLQ 

27 

0 

.0728865824, 

.0747231691,  .0761103619,  .0770398182,  .0775059480) 

GLO 

28 

X  «  A  *  G2  *  £2/ (2. Ml.  ♦  F NU2) ) 


40—NU0E~GAUSS'lT(ri3uA'0RA'rUR'E-F0'R'— INTEGRAL' AL ONG  FIBER 

TO  «  (2*8  -  RH01  /  2 .  t  T I  *  RHO  /  2.  i  FL  -0. 

00  10  N*1 1  20 .  . . . 

TEMPI  *  TO  -  T l *X4 ( N )  i  TEMP2  «  TO  ♦  T1*X4(N) 

CALL  TAUXZTTX  .“TEMPI,  ~ANST)  . . . 

CALL  TAUXZ1IX,  T  EMP2 ,  ANS2I 

FI  *  F  l  ’♦  AVflY)  '  *  "  (  ANS  I  "V  ANS2 1  . 

10  CONTINUE 

FI  «  Fl  '*  TV  .  . 

C  fo“ nOu E  CaTTs" S Tan  QTTa'oRaT  u R F“ INTEGRA l~‘a long  resin 

c 

TO'  *  -RH0”/”2.' . *V . tf  -  2  .  B  ♦  K  HO )  /  2.'  $  F2  *  0. 

r. 

DO  20  N=  1, V  . " . .  . . . 

TEMPI  «  TO  -  T 1  *  XUN) 

TEMP2  *  ’TO  ♦  Tt  *  XlTN) 

CALL  TAUX/2IX,  TEMPI,  ANSI) 

CALL  TAUXZ2('x,'YEMP2,~  ANS'2)~ .  . 

F2  »  F2  ♦  MINI  *  (ANSI  ♦  ANS2) 

20'  CONT  I  NUE  .  . . 

_ F2  ■  F_2  *__TJ _ 

C 

SHR  «  (FI  ♦  F 2 )  /  2. 

GCOMB"  *  SHR  *SGRTF(3.)  v  . 

GG2  *  GC0MB/G2 

p  r  i  n  r '{.  oo Vf" iV'ftT's  hr',"GC  om'b7‘  gg"2  . 

c  _ _ _ 

600  FORMAT'!  3H  fl ,  1  2X ,  E 1  5 . 4 ,  /  ,  3H0T2 ,  lTx ,  E 1 5 . 4 ,  /  SHOYbAR  ,Tox,'E15 .  V,7',“ 
*6H0GC0MB,9X,E15.4,/,5H0G/G2,9X,E IS. 4,/, 1H1 ) 

. . "RETURN  i  END  . 
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ooo  ‘  •  0  0-0  :  o  .000! 


SUBROUTINE  GRID  . . 

COMMJN/INPI/A,B,RHO,E_l,F»JUl,E2j_fNU2_  _  _ 

'  COMMON/  PO|"ntY/ XXTlVO  I  .VyiTSO  I  ,'JFL  AGfl  50T 
COMMON/LABL/  1N0EXI200) 

C 

DIMENSION  I Y l ( 7 ) i  1X1(6),  IY2I7),  1X2(6) 

OATA(IYl'«  0,  6,To,"YaV"  16,'  18,  20) 

_  DAT  A  (  I XI  »  1, _ 5_, _ 7 , _ 1 1 ,  IS) _ 

DAT  A (  I Y2  •  0,  A,  6,  d,  10,  14,  lb) 

0ATAIIX2  «  17,  13,  9,  S ,  3,  1)  _ 

C 

FXIX)  -A  -  SORTF ( RHOSO  -  IX-B)**2) 

F YI  X )  *'  B  -“SORTF  (  RHUSb’-  IX- A )  *«2  ) 

C _ _ 

RHOSO  «  RMO  •  RHO 
J  -  1 

OX'  ■  A/  19 :  $  DY"i  e/lY. 

'*•*  points'  Vn "Fiber . . . 

00  To- i~6 

IlY  -  lYlfK)  »  I2Y  -  IYI  (K  +  l  )  -  2 

I IX  VixiTkY . 1  lix"«  l 7  ~~ 

00  5  IY  •  11Y.I2Y.2  _i  FFY  -  IY 

Y  *  fc"  -  F  F Y”*'”dV . 

00~s~Tx~~lTxTlTx75  i  f  F  x*“  fx 

X  »  DX  *  F FX  _  _ 

. x'XljfV*  ■»  UlTfr Y  »'  JFL AGtJI  •*! 

J  -  J  ♦  1 

S  CONT I  NUt:  . 

10  CONTINUE _ 


**_* _ POINTS  IN  RESIN _ _ _ _ _ 

DO  20  K«  1 ,6 _ _ _ _ 

. IlY  V[  Y2<X)  '  'V  '  *  I2V  ■  7y2TkVi)"‘-  2 

_ I  _ I _ L£i_=_lMl><J _ 

DO  15  IY  -  1 1 Y,  I2Y,  2  i  FFY  «  IY 

Y  «  FFY  *0?  -  (J 


DO  15  IX  ■  IlXx  12Xj  }_  _  » _  FFX  »  IX 

X  «  OX  *  FFX 

_ PUL-  x _ I _ Y*U)  ■  V _ i _ AgiMLU  ■  i 

J  -  J  ♦  1 

15  continue 

20  CONTINUE 

'  ■'  po  i  NTS*"aLonT7nTe'aTa'ce'  . 


Y  ■  B  » 

XXI 109)  «  X  * 

Y  »  B-  2 .  *0Y  t 

XXI 110)  «  X  1 

x  ■  ux'  *  i 

_ _ L 

Y  -  B-  4. *DY  t 

xxi  li?i.  -  x  _ * 

Y  ■  B-  6 . *DY  $ 


X  ■  A-RHO 
YY(  109)  -  _Y 

X  ■  FXIY) 
YY ( 1101  *  Y 
**"  V  FYlx'l 
M.(UI)  •  1 _ 

X  ■  FXIY) 


YYI112)  «..Y 

X  ■  FXIY*) 


JllJLlAi-S-I 


t 

» 

1 

t 


.JFLAGI  109.)  .  »  3... 
JFLAG(llO)  «  3 
■  JEk*6t  1.UL1-J— j__ 
JfkACAM?)  -  \ 
JFk.*PAlAjJ—-— A _ 
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JFLAGI11A)  ■  3 


V  *  B-  d.*OY  t  X  ■  FXIY) 

XXI l IA)  «  X  i  YY( 11A)  »  Y  » 

X  «  3.*  OX  *  Y*v  F  Yt'Xl 


xxi 1 18)  *  x 

YYI  1161 

*  Y 

* 

JFLAGI 1151 

« 

3 

Y  *  B-10.*DY 

X 

*  FXIY) 

XXI 1 16)  *  X 

YYI 1 16) 

■  Y 

S 

JFLAGI 116) 

S 

3 

Y  »  H-12.*0Y 

X 

«  FXl  Y  1 

XXI  1  t  n  *  X 

YYI 1  17) 

«  Y 

t 

JFLAGI 1 17' 

X 

3 

x  *  5.*nx 

Y 

«  FYIX) 

XXI  1  Id)  •  X 

YYI 118) 

*  Y 

s 

JFLAGI 1  Id) 

* 

3 

Y  «  B-lA.*l)Y 

X 

*  FXIY) 

XXI 1 191  <  X 

YYI 119) 

*  Y 

A 

JFLAGI 1 19) 

= 

3 

x  *  7.*nx 

Y 

*  FYIX) 

XXI120)  *  X 

YYI 120) 

•  Y 

t 

JFLAGI 120) 

a 

3 

X  «  9.*DX 

Y 

«  FYIX) 

X  X  (  1  2  1  )  *  X 

YY I  121  ) 

*  Y 

* 

JFLAGI 121 ) 

a 

3 

Y  «  (3-16. *0Y 

X 

*  FXIY) 

XXI 122)  *  X 

YYI 122) 

a  Y 

% 

JFLAGI 122) 

a 

3 

X  *  11. *ox 

Y 

*  FY  1  X  ) 

XX  I  1  2  3  )  *  X 

YYI 123) 

*  Y 

% 

JFLAGI 123) 

a 

3 

X  *  l  3.  *0X 

Y 

«  F  Y  Ox  ) 

XXI  l  2  A  )  «  X 

YYI  1 2  A  ) 

a  Y 

% 

JFLAGI 12A) 

a 

3 

V  «  B-1M.*DY 

"X 

*  FXIY) 

XXI 125)  *  X 

YYI 128) 

«  Y 

l 

JFLAGI 125) 

a 

3 

X  •  15.*DX 

Y 

-  FYIX) 

XXI 126)  *  X 

YYI 1261 

■  Y 

% 

JFLAGI 126) 

a 

3 

X  «  1  7 *  0  X 

Y 

*  rr rn 

XXI 12  n  «  X 

YYI 1271 

*  Y 

t 

JFLAGI 127) 

c 

3 

X  «  A 

Y 

«  Q  -  RhO 

XXI 1281  »  X 

YYI  128) 

■  Y 

1 

JFLAGI 1  2d  I 

X 

3 

I NUt  X { X  ?H )  «  “500 

c 

no  "30  J«'l , 1*2> 

INOFX(J)  *  .J 

30  CONI INUt 

c 

J  •  129 

X  -  A  i  III  ■  0 

00  AO  lYVUY.lf  i  FFY  •  2*IY 

Y  «  B  -  FFY  *  OY 

x  x  ( j  ♦  i  y  )  *  *  x  V  yV(  j'+Ty  V  V  y" . i . -  •  - 

JFLAG(JMY)  *  1  t  _ 1N0EX1 J* I Y >  «  1000  ♦  IV 

AO  CONTINUF 
C 

JF  LAG (139)  •” J  F L A G (l AOV  *  2 
JFLAG(lAl)  •  A 
C 

RETURN  i  fcNO 


CARO  TOTAL  *~T22? 


BLANK  PAGE 


Unc  lassified 


Security  Climricition 


DOCUMENT  CONTROL  DATA  .RAD 

(Security  ot  tltlm.  body  of  abstract  and  Indaatna  annotation  must  bo  antatad  rfcwi  the  am  -W  ranart  la  alaaati 


I  OR1QINATINO  ACTIV  TV  {CtHpatOtO  author)  It#.  REPORT  SSCUMITV  CLASSIFICATION 

Whittaker  Corporation  |  Unclassified 

Research  h  Development/San  Diego 
San  Diego,  California 


9  RIRORJ  TITL* 

OBLIQUE  LOADING  OF  UNIDIRECTIONAL  FIBER  COMPOSITES; 
SHEAR  LOADING 


4.  ocscfiftivk  NOT*!  (T Yr*  •>  rtpotl  mti  Inehmln  *>»«; 

Final  Report  (March  1967  -  May  1968) 


t.  >utHO!{llfnnlMM,  mlddtm  Initial.  Ia> Inrnr' 

Juan  Haener 

e 

Alberto  Puppo 
Mipg-Yuan  Feng 


!  FtFOFT  OAT* 

January  1969 


M.  CONtHACT  an  CHANT  NO. 

Contract  DAAJ02-67-C-0035 

».  MOJIC  T  NO 

Tark  1  F  1 62204  A  1700- 


IT*.  TOTAL  NO.  OF  !AHI  Tk.  HO.  OF  NC  FI 


m.  oniaiNATom  ncfoft  nui/binui 


USAAVLABS  Technical  Report  68-81 


ft*,  other  RKftORT  NOO)  (Any  othar  ni mboro  Gat  may  bo  aaalgnad 
thta  roport) 

Whittaker  MJO  319  Final  Report 


10.  Dt  I  Tfll  BUT  ION  ITATIMINT 


This  document  has  been  approved  for  public  release  and  sale;  its  distribm;on 
is  unlimited. 


II.  IPOMIOIIINI  MILITANT  ACTIVITY 
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Stress  fields  in  the  components  of  a  multifiber,  unidirectional  composite  that 
are  caused  by  shear  loading  in  axial  directions  are  studied,  and  solutions  are 
presented.  Computations  for  various  volumetric  contents  and  material 
constants  are  described,  as  well  as  experiments  which  showed  good  corre¬ 
lation  with  the  theoretical  results  from  this  work.  Further,  it  is  shown  that 
certain  mathematical  solutions  normally  used  for  such  problems  are  not 
compatible  with  present  boundary  conditions,  and  that  the  problem  becomes 
two-dimensional . 
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